R.P. I. Technical Report MP-30 


A Progress Report for 
July 1, 1972 to December 31, 1972 

ANALYSIS AND DESIGN OF A CAPSULE 
LANDING SYSTEM AND SURFACE VEHICLE 
CONTROL SYSTEM FOR MARS EXPLORATION 

National Aeronautics and Space 
Administration 

Grant NGL 33-018-091 


Submitted by the Special Projects Committee 

D. K. Frederick 
P.K. Lashmet 
G.N. Sandor 
C.N. Shen 

E. J, Smith 
S.W. Yerazunis 


School of Engineering 
Rensselaer Polytechnic Institute 


ABSTRACT 


Investigation of problems related to the design and control of a 
mobile planetary vehicle to implement a systematic plan for the explora- 
tion of Mars has been undertaken. Problem areas receiving attention 
include: vehicle configuration, control, dynamics, systems and propul- 

sion; systems analysis; terrain modeling and path selection; and chemical 
analysis of specimens. The following specific tasks have been under 
study: vehicle model design, mathematical modeling of dynamic vehicle, 

experimental vehicle dynamics, obstacle negotiation, electromechanical 
controls, collapsibility and deployment, construction of a wheel tester, 
wheel analysis, payload design, system design optimization, effect of 
design assumptions, accessory optimal design, on-board computer, subsystem, 
laser range measurement, discrete obstacle detection, obstacle detection 
systems, terrain modeling, path selection system simulation and evaluation 
gas chromatograph/mass spectrometer system concepts, chromatograph model 
evaluation and improvement. 

These tasks are defined in detail and the progress which has been 
achieved during the period July 1, 1972 to December 31, 1973 is summarized 
Projections of work to be undertaken in certain areas during the period 
December 31, 1972 to June 30, 1973 are included. 
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Analysis and Design of a Capsule Landing System 
and Surface Vehicle Control System for Mars Exploration 


I. Introduction 

Current national goals in space exploration include a detailed study 
and examination of the planet Mars. The effectiveness of Mars exploration 
missions would be enhanced according to the extent to which the investigative 
devices which are landed are mobile, to the range of their mobility, and to 
the ability to control their motion. In order to achieve basic mission objec- 
tives, and beyond that, to maximize the return on the commitment of resources; 
to the mission, formidable technical problems must be resolved. The major 
factor contributing to these problems is the round trip communications time 
delay between martian and earth control stations which varies from a minimum 
of about 9 minutes to a maximum of 25 minutes. This time delay imposes strin- 
gent requirements on the vehicle, on its control and communication systems and 
on those systems included on board to make the scientific measurements, in 
terms of their ability to function autonomously. These systems must be able 
to operate with a high degree of reliability and must be capable of calling 
for earth control under appropriate circumstances. 

A number of important problems originating with these factors and 
relating directly the basic mission objectives of an unmanned exploration of 
Mars have been and are currently being investigated by a faculty- student pro- 
ject team at Rensselaer Polytechnic Institute with the support of NASA 
NGL-33-018-091. 

This progress report describes the tasks which have been undertaken 
and documents the progress which has been achieved in the interval July 1, 

1972 to December 31, 1972. 

II. Definition of Tasks 


The delay time in round trip communication between Mars and Earth 
gives rise to unique problems relevant to martian and/or other planetary 
explorations. All phases of the mission from landing the capsule in the 
neighborhood of a desired position to the systematic traversing of the sur- 
face and the attendant detection, measurement, and analytical operations 
must be consummated with a minimum of control and instruction by earth based 
units. The delay time requires that on board systems capable of making 
rational decisions be developed and that suitable precautions be taken against 
potential catastrophic failures. Four major task areas, which are in turn 
divided into appropriate sub-tasks, have been defined and are listed below. 

A. Vehicle Configuration, Control, Dynamics, Systems and Propulsion . 

The objectives of this task are to investigate problems related 
to the design of a roving vehicle for Mars exploration with 
respect to configuration; motion and attitude control; obstacle 
avoidance; control, information and power systems; and propulsion 
systems. In addition, the design concepts must accommodate the 
equipment and instruments required to automate the vehicle and 
to perform the scientific objectives of the mission. 
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B. General Systems Analysis . The objective of this task is to 
develop a framework within which decisions in design involving 
conflicting requirements can be made rationally and in the con- 
text of the whole system and mission. Relationships between 
alternative mission profiles and specifications and weight, 
energy and space allocation and management will be sought. 

C. Surface Navigation and Path Control . Once the capsule is landed 
and the roving vehicle is in an operational . state, it is neces- 
sary that the vehicle can be directed to proceed under remote 
control from the landing site to specified positions on the 
martian surface. This task is concerned with the problems of 
terrain modeling, path selection and navigation between the 
initial and terminal sites when major terrain features preclud- 
ing direct paths are to be anticipated. On board decision making 
capability must be designed to minimize earth control responsi- 
bility except in the most adverse circumstances. 

D. Chemical Analysis of Specimens . A major objective of martian 
surface exploration will be to obtain chemical, biochemical or 
biological information. Many experiments proposed for the 
mission require a general duty, gas chromotograph-mass spectro- 
meter for chemical analysis. The objective of this task is to 
generate fundamental data and concepts required to optimize this 
chemical analysis system. 

• III. Summary of Results 

Task A. Vehicle Configuration, Control, Dynamics, Systems and Propulsion 

This broad task has been divided into the following subtasks: vehicle 
model design, construction and evaluation /nathematical modeling of vehicle dyna- 
mics; experimental dynamics analysis; obstacle negotiation; electromechanical 
control systems; c'ollapsibility and deployment; wheel testing apparatus; and 
wheel analysis. A brief summary of the progress achieved in the prior inter- 
val in each of these areas follows below. 

A.l.a. Design, Construction and Evaluation of a 0.4 Scale Model 

A 0.4 scale model of the RPI-MRV is being constructed to serve 
a two-fold purpose. The first objective is to permit implementation 
of vehicle collapsibility and serve as a test bed for vehicle com- 

# ponents. The second objective is the comparison of computer simula- 
tion results with physical testing of the 0.4 scale vehicle. 

A.l.b. Vehicle Dynamics 

The work of this task is continuing in the development of a 
mathematical model of the RPI-MRV. The emphasis on analog simula- 
tion has been shifted to that of simulation on a digital computer. 
Preliminary results of the simulation have been used to determine 
spring rates for the 0.4 scale vehicle. In the longer range, this 
mathematical model as validated by experimental studies will provide 
the basis for optimizing the design. 


A.l.c, Stability and Obstacle Negotiation 

The programs developed through June 1972 for analysis of 
vehicle behavior over step, crevasse and slope obstacles were 
to be modified to include the adjustable parameter of payload 
height. In this effort, a force-moment approach was taken 
instead of the previously used geometric approach. This simpli- 
fied the analysis and will permit the development of a routine 
to evaluate stability and torque requirements for compound 
straight line obstacles and obstacle fields in a plane. 

A.l.d. Radio Control 

With the development of a 0.4 scale model to be used as a 
test bed for vehicle component study, the remote control of the 
vehicle is desired. This task has determined the mode and 
manner of signal transmission to link the vehicle w^th a remote 
controller. Hardware is being designed and built to implement 
the scheme developed. 

A. 2. Collapsibility and Deployment 

The vehicle has taken a new form with respect to its pay- 
load shape and location. This new form complements the collapsi- 
bility and deployment scheme which has been developed. The 
payload may now be located in the center of the capsule which 
means the placement of the center of gravity central to the 
capsule is only a problem of arranging payload components so the 
• payload eg is central to the payload. This condition is de- 

sirable for vehicle operation as well as meeting, capsule con- 
straints. 

A. 3. a. Wheel Tester and Grouser Design 

Three aspects of making the wheel tester operational and 
using it to test wheels have been completed. The soil which 
will be used to simulate the Martian surface, as it is under- 
stood, has been prepared. The instrumentation required to make 
the wheel tester operational has been developed. Grouser de- 
signs have been prepared as candidates . for testing. These 
designs incorporate variations in the basic parameters of wheel 
tread design. The data taken from experiments with these de- 
signs will be used to determine an optimum design for the 
conditions used in the testing. 

A. 3 . b . Wheel Analysis 

An effort was made to correlate past wheel design programs 
with existing models of the wheel. Correlation was not con- 
sidered acceptable, so a new approach was taken to the problem 
of wheel analysis for design purposes. Correlation of test 
data with the results of the new analysis is acceptable for 
small deflections. This method will be extended to include 
deflections which go beyond the linear range, an assumption 
which had been made to simplify the first stage of the analysis. 
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Task B. General Systems Analysis 

This major task is subdivided into three subtasks: identification 
of optimal designs and their sensitivity; time simulation of vehicle opera- 
tion; and specification of an integrated Mars Roving Vehicle controller. 

B.l. Identification of Optimal Designs and Their Sensitivity 

The system design optimization problem has been solved for 
the case of a 4-wheeled rover communicating directly with earth. 

. The solution specifies values of the parameters of the system 
design, e.g., weight and power allocations to the vehicle sub- 
systems, and subsystem performance levels (communications data 
rate, vehicle velocity, science time). In addition, two other 
cases (4-wheeled rover communicating via a Mars orbiter, 6- 
wheeled direct communicating rover) will be solved by July 1973. 

A method has been developed which identifies the optimal 
readjustments of design parameters to a forced perturbation in 
any one parameter of the design. This procedure will be used 
to investigate the sensitivity of the optimal designs identified 
above, and promises to be of use in dealing with the effects of 
any design variation. 

B.2. Time Simulation of Vehicle Operation 

Investigation of the operation of the vehicle system in 
real or scaled time was initiated during the reporting period. 

The purpose of this subtask is to analyze the operation of the 
designs identified in (1) with differing mission profiles and 
stochastic inputs, (terrain slopes, emergencies,...). The ef- 
fects of an increased amount of earth control can also be 
determined. A computer simulation in the Monte-Carlo mode will 
be prepared by May 1973. The simulation will "track" the vehicle 
status during its on-surface operations. 

B.3. Specification of an Integrated MRV Controller 

The problem of onboard computer control of the many vehicle 
operations is being investigated. Specifically, the demands to 
be placed upon the onboard computer by the vehicle subsystems and 
queueing and system management operations by the computer will 
be specified and simulated on a GE-PAC/30 which has been pro- 
grammed to have the capability of "appearing" to be any of a wide 
range of computational machines. A working simulation which 
includes a section which simulates the MRV as seen by the com- 
puter will be available by July 1973. 

Task C. Navigation, Terrain Modeling and Path Selection 

This task deals with problems of terrain sensing, terrain modeling, 
obstacle detection and path selection for an autonomous roving vehicle. 
Specific tasks receiving active effort include: the design of a laser range- 

finder for terrain modeling; terrain sensor instrumentation for a laser scan 
system; stochastic estimation of gradients; stereo and non-stereo obstacle 


detection systems; end path selection systems evaluation. 

• C.l.a* Design of a Laser Rangefinder for Terrain Modeling 

Various methods by which a laser may be used for rangefind- 
ing purposes in the range of 3-30 meters have been analyzed and 
evaluated. Optical focussing, timed puise and phase difference 
techniques were compared. It has been concluded that a hybrid 
system combining the timed pulse and phase difference concepts 
is most promising. The construction of a system to be used 
under laboratory conditions to determine reflectance of typical 
materials and their accuracy is currently underway. 

C.l.b. Terrain Sensor Instrumentation - Laser Scan Systems 

. The problem of obtaining scan information in the 3-30 meter 
range has been considered in detail. Emphasis has been directed 
to electro-optical methods which would permit scanning at a 
sufficiently high rate to allow for vehicle motions to be neglec- 
ted. Acusto-optic beam benders, piezoelectric deflectors and 
digital light deflectors were considered. It was concluded that 
a system based on digital light deflectors was the most feasible. 
A detailed scanning system has been outlined. 

C.l.c. Stochastic Estimation of Gradients 

An analysis of the effect of measurement errors, i.e., 
range and elevation and azimuthal angles, on the gradient of 
terrain in the 3-30 meter range has been undertaken. It is found 
that measurement errors originating with vehicle motion, (i.e. 
roll, pitch and heave) result in very large errors in estimating 
the in-path and cross-path slopes and the gradient. It is con- 
cluded that a rapid scan system should be used to eliminate these 
errors. Assuming a rapid scan system obviating vehicle motion 
effects, the relationships between the standard deviation of 
gradient as a- function of range and of the standard deviation 
range and elevation and azimuthal error for, a horizontal flat 
plain have been determined. The standard deviation of gradient 
varies from less than 1° to about 8° depending on the measurement 
errors. The effect of an inclined terrain has also been deter- 
mined. 

As this effort more completely relates the effect of measure- 
ment errors on the reliability of terrain model interpretation, 
specification for a terrain sensing system can be defined more 
rationally. 

C.l.d. Stereo and Non- Stereo Obstacle Detection Systems 

Stereo and non-stereo obstacle detection concepts have been 
evaluated with regard to their ability to identify the height 
and range of discrete obstacles. The three systems studied were: 
(a) stereo angles, (b) stereo range, and (c) non-stereo (range- 
angle). It is concluded that the stereo angles method which 


involves only the measurement of two angles is feasible only if 
the range measurement error cannot be reduced below 10 cm. and if 
the range is less than 12 meters. The stereo range method can only 
be effective if range errors can be reduced significantly. The 
non-stereo method is almost as good as the stereo angle method 
at short range and is much better at ranges larger than 12 meters. 

C. 2. Path Selection System Simulation and Evaluation 

Development of the computer program for realistically simula- 
ting the performance of various path selection systems has con- 
tinued. A variety of improvements have been incorporated including 
the vehicle's dimensions and attitude, both vehicle and vertical- 
fixed sensors, improved range measurement resolution, and improved 
output capabilities. 

A limited comparison of three related path selection systems 
which use single and split beam laser sensors has been carried out. 
Currently work is being devoted to refinement of the program and 
the evaluation criteria, in addition to a comparison of vehicle- 
and vertical-fixed sensor configurations. 

Task D. Chemical Analysis of Specimens 

This task is concerned with developing the fundamental concepts 
which will be required to optimize a gas chromatograph-mass spectrometer 
analysis system. Topics receiving active effort include mass spectrometer 
characteristics, carrier gas generation and removal, multicomponent chroma 
tography, and chromatograph model improvement. 

D. l Gas Chromatograph/Mass Spectrometer Systems Concepts 

D.l.a. Mass Spectrometer System Characteristics 

The mass spectrometer studies are intended to characterize 
the major properties of the instrument, the resulting relations 
to be suitable for investigation of system power, weight, and 
volume requirements. Relations involving mass resolution for 
single and double focussing instruments have been developed. 
Preliminary studies showed the theoretical analysis of instrument 
sensitivity to be beyond the scope of this project, and empirical 
relations are being sought. Investigation of the effect of major 
design variables upon power requirements has commenced with em- 
phasis being directed on the ionization and accelerating chambers. 

D.l.b. Carrier Gas Generation and Removal 


Preliminary analysis of high-pressure gas storage as a source 
of carrier gas has been completed. Use of filament reinforced 
construction will provide storage capacity for one year's opera- 
tion of the chromatograph within a volume of about one or two 
liters and a weight of about 2 kgm. The concept appears competi- 
tive with a previously proposed electrolytic gas generator. En- 
gineering analysis of carrier gas removal by diffusion through a 
palladium alloy tubing has commenced. 
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D.2 Chromatographic Systems Analysis 

D.2.a. Multicomponent Chromatography 

Prior results showing limitations on superposition of pure 
component data to predict multicomponent chromatograms have been 
verified. Additional multicomponent data are being obtained on 
porous and nonporous adsorbents. Improvements in the chromato- 
graphic test facility are being undertaken to permit more auto- 
matic processing of the data. This involves magnetic tape 
recording of the composition detector signals and preliminary 
data processing by minicomputer. Experimental characterization 
of the adsorbents (pore size distribution, etc.) has begun. 

D.2.b. Chromatograph Model Improvement 

Intraparticle diffusion and rates of adsorption have been 
incorporated into the mathematical model of the chromatograph 
to improve representation of the experimental data. Moment 
analysis of the resulting system of equations shows the model 
to be capable of representing more adequately the spreading of 
experimental data. It appears an analytical solution of the 
model is not realizable, and numerical solution methods are be- 
ing investigated. 

IV. Detailed Summaries of Progress 

„ Task A. Vehicle Configuration, Control, Dynamics, Systems and Propulsion 

This broad task is concerned with all aspects of the design of a 
roving vehicle for Mars exploration. Included within this major effort are 
the following subtasks: vehicle model design, mathematical model of the 

vehicle dynamics, experimental dynamic analysis, obstacle negotiation, 
electromechanical control systems, collapsibility and deployment, wheel test- 
ing apparatus, wheel analysis, and payload design. 

Progress made in these areas during the past interval is summarized 
in the sections which follow. 

A. 1. Vehicle Design, Construction and Evaluation 

Student Group Leader: J. Almstead 

Faculty Advisor: Prof. G. N. Sandor 

Areas of major concern to the Vehicle Group are vehicle 
control, component design, dynamic analysis, and vehicle evalua- 
tion with regard to mission parameters. Vehicle control pertains 
to the basic functions of the MRV: steering, drive speed, pay- 

load attitude, and vehicle configuration. The problems associated 
with the vehicle control include sensing the states, acting upon 
and coordinating the functions. Component design pertains to the 
development of what might be considered subsystems of the vehicle 
structure. The major emphasis in this area has been placed upon 
payload and wheel design. Dynamic analysis has, and will include 
both analytical and experimental approaches to the determination 
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.of the dynamic behavior of the MRV. Vehicle evaluation pertains 
to determining the performance factors of the MRV and comparing 
them with the performance criteria of the mission objectives. 

The Vehicle Group has developed what is considered a third 
generation design of the RPI-MRV. The current design includes 
certain payload and collapsibility aspects which the previous 
models did not include. The new features have changed the shape 
of the vehicle slightly, but retained the original concepts of a 
dragster vehicle with wagon steering. The vehicle design has 
reached a stage where evaluation with regard to the mission goals 
and the capabilities of other proposed vehicles is necessary. As 
the evaluation is conducted, areas requiring additional redesign 
will become apparent, and those adjustments will be made. 

. The construction of a 0.4 scale model of the vehicle will 
permit physical testing and evaluation of the design capabilities 
and limitations. The model, presently near completion, will be 
ready for testing during the spring semester. The analysis of the 
design and testing of the model will be directed toward establish- 
ing performance factors in the areas of dynamic behavior, stability 
and obstacle negotiation, maneuverability, payload fraction, flexi- 
bility, reliability, and capsule fit. The plans for remote and 
computer control of the model will provide the vehicle the capa- 
bility of being a test bed for the analysis of vehicle control 
concepts. 

Projects planned for the spring semester are in the areas of 
dynamic analysis, wheel development and testing, vehicle control 
and design evaluation. 

A. l.a. Construction of a 0,4 Scale Model - William Embleton 

Faculty Advisor: Prof. G. N. Sandor 

The objective of this task is to design and construct a 0.4 
model vehicle- structure which would support a scaled load comparable 
to the full size vehicle, which would be collapsible and maneuver- 
able without an umbilical cord. According to Ref. 1, the estimated 
weight of a full size vehicle is projected as 1100 pounds. A com- 
parable scaled weight for the 0.4 model was calculated to be 70 pounds. 

An area moment analysis was performed on the rigid members of 
the vehicle to identify a configuration of each part which would 
allow them to be light and strong. Assuming a maximum static deflec- 
tion for the torsion bar and for each of the structural members, 
appropriate dimensions of each part were estimated. 

While the linear wheel model was programmed, a non- powered 0.2 
model was built to finalize the vehicle design to be incorporated 
in the 0.4 model. The 0.2 vehicle model implemented the collapsibi- 
lity scheme. It also utilized the double torsion bar which allows 
more room in the payload by putting the levelizing and attitude motors 
on the main frame. 

Construction is now in progress on the 0.4 model and it is 


planned that the construction of the frame with the motors will 
be completed by the first week of February 1973. When completely 
finished the 0.4 vehicle model will be collapsible, radio con- 
trolled, and self-powered. During the next period, . dynamic test- 
ing of the 0.4 model will be initiated. 

A.l.b. Vehicle Dynamics - Mark C. Rodamaker 

Faculty Advisor: Prof. G. N. Sandor 

The objective of this task is the mathematical modeling and 
computer simulation of the vehicle to predict vehicle response 
to perturbations. This modeling will permit estimation of the 
design parameters to insure a dynamic behavior of the vehicle 
which is compatible with the mission objectives. The modeling 
is being done assuming a bicycle configuration (Figure 1), and 
will be extended, in the future, to a three dimensional model. 

Equations of motion for the seven-degree-of-freedom planar 
model have been derived using a complementary virtual work method 
. which allows the inputs to be displacements rather than forces as 
required in classical Newtonian methods. An analog simulation 
was attempted, but due to problems with the analog machines avail- 
able, it was decided to discontinue the analog simulation and 
shift the emphasis to a digital simulation. The equations have 
been programmed, but at this time the program is not running 
correctly. 

Other work the past six months has included determining de- 
sirable spring constants for the wheels and torsion bar of the 
0.4 scale model. These values were based mainly on analysis of 
test data accumulated over the past two years on the previous 0.184 
scale models. 

Emphasis for the next period will be on the digital simulation 
of the s.even- degree- of- f reedom model and investigation of the de-. 
velopment of a three-dimensional simulation. These results should 
provide a rational basis for improving the design of the 0.4 scale 
model. 

A. l.c. Stability and Obstacle Negotiation - Larry Mabius 

Faculty Advisor: Prof. G. N. Sandor 

Previous work in the area of stability and obstacle negotia- 
tion concentrated on determining the rear-wheel torque as a 
function of displacement as the obstacle was negotiated. Earlier 
studies assumed a fixed vehicle configuration and used trigonometric 
analysis to calculate the torque values. 

Because the vehicle is capable to adjusting its payload height, 
it was felt that the vehicle configuration should be made an adjust- 
able parameter. To facilitate this, a force-moment analysis was 
used in calculating the torques required. With the new approach, 
the vehicle configuration together with the three angles indicated 
in Figure 2, suffice for computing torque requirements. One angle 



Figure 1. Mathematical Model of MRV 



Is the pitch of the vehicle and the other two reflect the slopes 
of the terrain on which the wheels are located. Also included 
was an indication when the center of gravity fell outside the 
wheel contact points, causing instability. 

Figure 3 illustrates the relation among the dimensional 
parameters as the payload height is varied. The parameters of 
interest are the height of the center of gravity, the horizontal 
distance from the center of gravity to the rear wheel, and the 
overall vehicle length U,P and Lt respectively. Variations of 
these parameters have a simultaneous effect on the required torque. 
Analytical results indicate that, when negotiating step and crevasse 
obstacles, peaks in the torque demands occur when the front wheel 
starts over the obstacle and again when the rear wheel starts over 
the obstacle. In the case of slope obstacles, the analytical re- 
sults, Figure 4, indicate that sharp peaks in the torque require- 
ments do not occur. 

Increasing the height of the center of gravity, which effec- 
tively shortens the rear strut length, causes a decrease in the 
torque requirement as the front wheel contacts the obstacle, and 
increases the torque demanded as the rear wheel contacts the 
obstacle. This may be explained by considering the vehicle as a 
class two lever. The wheel going up at the step applies the 
e ff°rt while the other acts as a fulcrum and the weight acts as 
the load. With the front wheel applying the effort, moving the 
center of gravity up and therefore rearward increases the ratio 
of the moment arms thereby decreasing the force applied. On the 
other hand, the same action with the rear wheel applying the 
effort, decreases the ratio of the moment arms thereby increasing 
the force. Since the force applied is generated by the torque on 
the rear wheels, 'the predescribed results take place. Therefore 
it seems the shifting of the center of gravity back and forth has 
a greater effect than shifting it up and down. This can be veri- 
fied by observing the change in torque at U equal zero where the 
horizontal change in the eg is zero with respect to U and the 
height change is U. The change in torque is zero in all cases 
verifying the above effect. 

Present work with the vehicle is to generalize the obstacle 
field. This field would consist of a series of straight lines 
connected together. The planar vehicle would traverse it in 
discrete steps. Then at each point the torque would be calculated 
along with an instability factor and a check for the hang-up. The 
instability factor would indicate where, along a line going be- 
tween the two contact points on the wheels, the eg would be pro- 
jected down on. In percentage terms, if this point were at the 
center of the line, it would be zero percent. Likewise, if it were 
at the end of the line (the eg would be right above the contact 
point), it would indicate 100% instability. The hang-up condition 
would test for contact between the ground and the terrain between 
the two wheels. This system, extended to three dimensions, will 
demonstrate the optimum vehicle height for negotiating any given 
terrain. 

The generalization of the obstacle field in the plane should 

















be completed by March 15, 1973. Extension of the obstacle field 
to three dimensions should be concluded by May 15, 1973. 

A. l.d. Radio Control - Mark Rinaldi 

Faculty Advisor: Prof. G. N. Sandor 

The objective of this task is the design and construction of 
a radio-control system for the 0.4 model of the RPI-MRV. Design 
of the overall system is divided into two parts: 

1. - The design and construction of the communications 

equipment that will transmit information to the 
vehicle. 

2. The design of the accompanying vehicle control 
systems. 

Attention was focussed primarily on the first of these two 
goals. Analog methods were rejected for signal transmission be- 
. cause of the difficulty in interfacing an analog system with the 
computer. The benefits of a digital transmission are primarily 
its large data handling capabilities and the relative ease with 
which it can be interconnected with the computer. 

There are several ways in which digital . signals can be trans- 
mitted in a serial fashion. The two which were considered are 
shown in Figure 5. In all serial systems, the transmitter and re- 
ceiver must be synchronized so that the receiver can correctly 
interpret the incoming data. In the bipolar return-to-zero signal, 
there is a zero voltage, or "dead space" between each data bit. The 
appearance of a dead space indicates to the receiver that one bit 
has been completed and that another will soon be transmitted. The 
advantage of simplicity is readily apparent in such a system. Un- 
fortunately such a bipolar signal cannot be handled by transmitters 
such as walkie-talkies. 

Two types of transmitters were considered. The analog radio- 
control type, which is similar to those used for flying model 
airplanes, did not provide the desired system flexibility. The 
logical alternative was the walkie-talkie, which has the advantage 
of a potential for two-way, although not simultaneous, transmission. 

To complement the walkie-talkie, the unipolar synchronius 
type signal format was chosen. With this signal, synchronization 
can be provided in two ways : 

1. The simultaneous transmission of a synchronizing 
signal (also referred to as a clock) along with 
the data. 

2. If the data rate is slow, a start pulse can be 
transmitted prior to the command, signaling the 
receiver that data is about to be transmitted. 

Such a pulse would activate and synchronize a 
clock generator at the receiver. 
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The first case mentioned above would require' two trans- 
mitting channels or multiplexing of the data line and the 
clock. Because of this difficulty the second method was 
chosen. 

The next decision was the word length. For the actual 
data, a word of only four bits yielding 15 possible commands 
was deemed necessary. However in addition to the four data 
bits, one parity bit, one flag bit and two spare bits (for 
possible expansion of the system) were added. This brings the 
total bit count per command to eight. A sample command is 
shown in Figure 6. 

The parity bit is used as a possible means of error detec- 
tion. The first seven bits of the command are examined at the 
transmitter for even parity (an even number of "l"s). If even 
parity is detected then the parity bit is a "0"; if odd parity 
is detected the parity bit is a "1". The addition of this bit 
causes the now eight bit command to have even parity. When the 
command is received at the vehicle, it is examined for even 
parity. If this condition is met, the word is accepted as a 
valid command, if not (if odd parity exists) then the command 
is rejected, and a signal is sent back to re- transmit. Obviously 
an even number of errors will go undetected at the receiver. 

Work has been concentrated on designing the Command Encoder. 
A schematic of this encoder is shown in Figure 7. As of this 
date, no concrete design work has been done on the various blocks 
at the receiving end. The most important unit here is the 
Command Decoder, which translates all binary commands into con- 
trol signals for the various vehicle subsystems. The Command 
Decoder will operate in two modes as listed below: 

1. Real Time Command "Mode (RTC) : In this mode 

commands that are sent to the vehicle are decoded 
and executed immediately by the Command Decoder. 

This work is now in progress. 

2. Stored Command Mode: Here commands that are 

transmitted are not executed immediately. Instead 
the Command Decoder acts as a processor, routing 
the information to the Command Programmer (merely 

a memory unit). At some pre-determined time, these 
commands are fed to the Command Decoder which then 
decodes them and executes the necessary functions. 

This mode would simulate those times when an actual 
vehicle would be out of radio contact with Earth. 

(This mode will be implemented in the future). 

A brief description of the RTC mode is as follows: 

When the vehicle received the START (ST) pulse, it is told 
that data should follow at approximately 5 ms (microseconds) later 
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This gives the receiver time to enable the inputs to the 
•Decoder storage register and also to synchronize the onboard 
clock with that of the transmitter. After 5 ms the data 
. starts coming in at a rate of 100 bits/sec. The first bit 
is the FLAG bit, and after 8 clock pulses the data completely 
shifts into the Decoder, and the FLAG is detected in the last 
register space (shown on the preceding page) the clock and 
data inputs are inhibited (to prevent further shifting of the 
data).. The parity detector then samples the word, and accepts 
it if the parity is even. If odd parity is detected the com- 
mand is rejected and a RETRANSMIT command is sent back to the 
transmitter (assuming two-way communication). Once the word 
has been accepted, bits 2-5 are decoded, executing the command. 

Of the 15 commands available, at least one will be reserved for 
addressing the Command Programmer. 

The project objectives for the reat of the year are: 

1. Construction of the proposed Command Encoder, 

March 15, 1973. 

2. Design and construction of the Command Decoder, 

March 15, 1973. 

3. Linking of the above units by radio communication, 

April 15,1973. 

4. Interfacing the Command Decoder with the actual 
vehicle sub- systems, May 15, 1973. 

A. 2. Collapsibility and Deployment - J. Almstead . 

Faculty Advisor: Prof. G. N. Sandor 

In the design of the RPI-MRV it is necessary that the aspects 
of collapsibility and deployment be incorporated and analyzed as a 
' functional part of the vehicle to insure a study of the complete 
vehicle system. The major objective of coliapsibility-deployment 
is to permit the vehicle to fit within the Viking capsule. The 
secondary objectives are to utilize the coliapsibility-deployment 
functions to further the flexibility of the operating vehicle, and 
to provide for collapsibility and deployment with reliability and 
minimum weight. 

A coliapsibility-deployment scheme was developed which pro- 
vides for the vehicle to fit within the Viking capsule, Ref. 2. 

The scheme requires that the vehicle be "folded" in three areas. 
Figure 8 indicates these areas of concern to be: 1) the front 
frame pivot, 2) the folding and lengthening of the rear struts, 
and 3) the required flexibility of the front axle-vehicle frame 
connection. 

The front frame pivot is to be used for coliapsibility- 
deployment only. Once the vehicle has landed on the Martian surface 
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and has been deployed to its operational configuration, the 
front frame pivot will be fixed. Collars placed over the 
pivoting fore-section will attach it rigidly to the rest of 
the frame. 

. The rear strut pivot and lengthening will be accomplished 
by a four bar linkage. The linkage was developed through the 
use of the complex-number technique of linkage synthesis for 
: three precision points and specified fixed pivot positions. A 
new computer program was written using this technique and used 
interactively on a graphic time sharing terminal to synthesize 
and analyze the four-bar mechanism that can be used to deploy 
the rear wheels. Once the rear wheels have been deployed, they 
will remain in the operational configuration for the rest of the 
vehicle life. 

The front axle— frame connection has been designed to permit 
collapsibility-deployment, minimize weight, and provide for 
vehicle symmetry. The front axle and wheels may be rotated at 
the front axle-frame connection about an axis parallel to the 
torsion bar. This maneuver is accomplished by the same motor 
used to steer the vehicle, and permits vehicle symmetry about 
the plane of the vehicle frame. This allows the vehicle to per- 
form an "emergency maneuver" without moving the vehicle body, 
swing the front wheels over the vehicle to the rear, and vice- 
versa. Thus reversed, the vehicle can drive away from a condition 
from which it could not have backed away. 

These three mechanisms required for collapsibility-deployment 
are being implemented and tested on both the 0.2 scale model and 
0.4 scale model. 

Work during the second .semester will be directed toward mak- 
. i- n g operational and utilizing the Integrated Mechanisms Program 
(IMP) from Dr. J. J. Uicker of the University of Wisconsin, Ref. 3, 
and which is to be modified for our purpose. IMP is designed to 
analyze two or three dimensional mechanisms having single or multi- 
ple degrees of freedom. Kinematic, static or dynamic analyses are 
performed resulting in the calculation of positions, velocities, 
accelerations, forces, natural frequencies and damping ratios for 
a given mechanism. It is hoped that the vehicle "mechanism" and 
its component mechanisms may be analyzed with IMP. Once operational 
this tool, it is expected, will aid the design as well as the analy- 
sis' aspects of the mechanical design associated with the vehicle. 

IMP is not complete. It will require writing several routines 
to adapt it to the RPI computer system. The program as received is 
not operational and will require debugging of some of its sub- 
routines. An exact timetable is impossible to establish because 
the magnitude of the work is yet unknown. Much of the work depends 
on how similar or dissimilar the computer systems are. 

3. a. Wheel Tester and Grouser Design - J. Ozimek and P. Marino 

Faculty Advisor: Prof. G. N. Sandor 

The purpose of this project is to determine the performance 
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aspects of wheels for a martian roving vehicle. The program was 
separated into six different segments, each of which had a 
specific goal. The first three segments were expected to be 
completed by January 1973 and the remaining three are expected 
to be finished by June 1973, The progress of the first three, 
i.e. sand preparation, instrumentation, and Grouser design, as 
well as an overall summary is as follows. 

SAND PREPARATION: It was desirable to simulate the soil 

conditions the wheels would encounter over a year's span on Mars. 

Of course, this meant that all possible combinations of soil con- 
could exist, from fine loose sand to hard bedrock. From 
a traction standpoint, loose sand represents the worst of possible 
conditions. For this condition , sand between two and sixty microns 
was obtained from a local creek bed. The sand was dried to zero 
percent moisture content. The moisture can be varied by the addi- 
tion of water in a controlled manner. The density of soil condition 
can be varied by compaction as well as the addition of rocks. 

There are few published results against which our findings can 
be compared. However, as long as soil conditions are kept consistent, 
valid differences can be determined between various tests and relative 
difference can be evaluated. 

The environmental conditions at the test site are not optimal. 
Since most testing will take place during the colder months, heating 
will tend to evaporate the moisture content of the soil. This 
necessitates continual testing of the soil in order to maintain 
acceptable levels of accuracy. The test conditions will be dupli- 
cated, if necessary, by comparison with original recorded levels of 
parameters. Moisture content can be measured by comparison of wet 
and dry weights of soil. This involves a simple process of eva- 
porating moisture content by means of a dryer. 

INSTRUMENTATION: There were three different types of instru- 

ments needed for measurement. One was needed to measure speed, 
one to measure draw bar pull and one to measure delivered torque. 

The speed of both wheel and test cart will be measured electrically 
and recorded on a strip chart recorder. A device for measuring 
speed has been built and will give readings that are recordable on 
a strip chart. There will not be a great variation in speed and 
this device will be accurate enough within the range desired. 

The measurement of draw bar pull can also be accomplished 
electrically and a device for doing this has been made. The re- 
sults will be fed into a strip chart recorder with speed being 
recorded on a second channel. 

There were problems with measuring the torque delivered to 
the wheel. Originally it was planned to measure the electrical 
power delivered to the motor driving the wheel. If the wheel is 
always traveling slower than the car (negative slip) this is an 
acceptable situation. However, if the wheel travels faster than 
the cart (positive slip) then power is being delivered to the motor 


by the wheel. In this situation it is necessary to measure 
mechanical power directly, since the motor (or generator) current 
is no longer readily correlated with the power input. 

The system chosen includes the use of four mechanical trans- 
ducers which will yield constant outputs. Two solid state tacho- 
meters will record cart speed and rotational speed of the wheel 
under test. The tachometers were designed to facilitate record- 
ing of data during testing (see Figure 9). Two strain gauge 
systems will yield draw bar pull on the cart and applied torque 
to the wheel. One objective of this system is to enable all four 
parameters to be observed and recorded continually. 

Several values will be kept constant for each run. These 
include cant angle of the wheel, load on the wheel, climbing angle, 
and sand conditions. Using these values and the four values being 
plotted, the characteristics of the wheel can be determined. The 
characteristics of most interest are the tractive coefficient and 
drive efficiency. Both of these values are typically functions of 
slip speed and vehicle speed, slip speed having the major effect. 
The tractive coefficient is the ratio of effective frictional force 
of the wheel on the sand to the normal load on the wheel. The 
tractive coefficient improves with more severe Grouser design. The 
efficiency, however, is highest with less severe Grousers as a rule. 
Therefore, the final Grouser design will involve a trade-off with 
respect to maximum power available from the vehicle and required 
performance characteristics. 

Calibration of the instruments will be necessary before and 
during testing. This is not difficult because all conditions be- 
fore actual testing are known or can be easily determined. The 
instruments can be kept within the same scale range and remain 
linear within the range of testing results. 

Consistency of performance will have to be determined from 
the test results. The number of test runs necessary to obtain re- 
liable data is important. Without any advance knowledge of how the 
device will perform, the number of individual tests remains unknown. 
When testing begins a large sampling will be taken and compared 
®8®-i^st itself to determine deviation. This should give some in- 
dication of how many tests for each condition is necessary and also 
determine the number of different conditions that can be tested 
within the allotted time frame. 

GROUSER DESIGNi The design of Grousers is a conceptual engin- 
eering problem. Few actual formulas exist for design of Grousers, 
most design has been accomplished by use of empirical formulas or 
off the shelf designs. These formulas give only qualitative ideas 
of how Grousers will perform under various conditions. Three 
separate designs have been suggested and will be compared against 
each other. Trade-offs can be made among the three in order to 
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gain an optimally functioning Grouser. 

The actual material used for the implements will be old tire 
rubber, because it is easily cut and shaped to desired forms. It 
will be glued on the rim in the present application but may be 
bolted on if necessary. 

Selection of the best design will be made by comparison of 
data and observations from level and slope tests. Of course, in 
this application only three designs are being tested and compared. 
Ideally, tests should be conducted over a larger variety of 
Grouser designs as well as trade-offs among them. 

FUTURE WORK: Three things remain to be accomplished: slope 

testing, level testing and lateral stiffness tests. It is ex- 
pected that level testing and slope testing will take up one 
month each. A lateral stiffness test is necessary for three- 
dimensional mathematical modeling. This will be the last part of 
actual testing. A written summary of results and conclusions will 
complete the project. 

A.3.b. Wheel Analysis - W. Embleton 

Faculty Advisor: Prof. G. N. Sandor 

The task of developing reliable methods for analyzing wheels 
for the Mars mission has been pursued for several years and a 
number of computer programs have been developed. However, none 
of the programs correlate effectively the experimental data for 
the 0.184 scale, 0.4 scale or full-size models of the RPI toroidal 
wheel. The objective of this task is to obtain an effective 
analysis method which can represent accurately the behavior of 
the wheel. An effective method will provide a rational basis for 
the design of wheels for defined specifications. 

A linearized math model of the wheel was made with the help • 
of Dr. C. A. Broniarek using methods outlined in Ref. 4. Pro- 
grammed on the computer, the math model gave the deflection of the 
wheel due to a static load applied at the hub of the wheel. In 
this math model, the inner hoops which are attached to the hub and 
to a rigid rim are linearized. A load applied at the hub causes 
the springs to deflect yielding a system of forces along the rigid 
rim. Assuming a flexible rim and applying the system of forces a 
little at a time to the rim causes the rim to deform. After each 
deformation of the rim is accomplished, the deformed rim is regard- 
ed' as rigid with the new deformation of the springs causing a new 
system of forces. When the sum of the forces of the springs caused 
by the deformed rim and the initial load is zero, an equilibrium is 
found for the wheel and a deflection of the wheel is found for a 
given load. 

The results from the computer program compared satisfactorily 
to the experimental results previously observed. The graph of the 
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full scale model wheel showing the comparison of the computer 
results and the experimental results is shown in Figure 10. An 
effort is now being made to provide for Grouser design in the 
linearized math model program and to design a non-linear math 
model to be programmed which would give dynamic responses of 
the wheel due to external loads which might be applied. This 
math model will divide external loads into incremented parts and 
apply them to a linearized math model which will have material 
properties and spring coefficients changed after each deformation 
so that each incremented load will use the properties at the pre- 
vious deflection to get an approximation for a deflection at the 
larger load thus giving a deflection curve for loads greater than 
those permitted by the linearized analysis. 

During the next period a non-linear math model of the wheel 
will be designed and programmed. This will give a better design 
tool for the wheel giving fairly accurate displacement-versus- 
load graphs for the dynamic modeling of the vehicle. 

Task B. Systems Analysis 

The objective of this task is to develop a framework within which 
vehicle design decisions involving conflicting requirements can be made. 

Work previous to this report period established the method of optimal system 
design for the Mars roving vehicle to consist of; developing a static system 
model, measuring the system performance as a function of model variables, 
and imbedding the optimization problem in a nonlinear programming (NLP) 
format and so locating the optimal design. In addition, study of the computer 
requirements of such a system were initiated. 

The progress of the work during this period will be summarized as 

follows : 

B.l. Identification of optimal designs for several possible 

system configurations, namely 

a. 4-wheeled vehicle, direct Mars-Earth communication 

b. 4-wheeled vehicle, communication via a Mars orbiter 

c. 6-wheeled vehicle, direct communication 

B.2. Analysis of the performance of the optimal designs obtained 

in B.l by Monte-Carlo scaled-time simulation. 

B.3. Specification of an integrated Mars roving vehicle (MRV) 

controller. 


B. 1. Identification of Optimal Designs 

Design of a vehicle intended to semi- autonomously rove the 
Mars surface while gathering scientific data, involves the in- 
corporation of the ability to perform eight rather distinct 
functions. These functional divisions, or subsystems are: 
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Science — collecting data and samples and performing 
experiments 

Communications - transmitting science and engineering 
(vehicle status) data to earth and receiving earth 
commands 

Power - generation and storage of energy required by 
vehicle subsystems 

Thermal control - maintaining acceptable internal tem- 
peratures onboard 

Navigation - locating the vehicle in Mars- centered co- 
ordinates 

Obstacle avoidance - identifying terrain hazards and 
choosing a safe path for vehicle travel 

Computation and data-handling - performing onboard 
computational operations 

Vehicle structure - frame, motors, and suspension 

The requirement that the vehicle system be deliverable to the 
Mars surface imposes weight, volume, and size limitations. The de- 
Problem then becomes an allocation problem. The operation of 
each subsystem will improve with increased power, weight, and 
volume allotted to it. The technique of this approach is to use a 
system model to measure these improvements and arrive at optimum 
allocations on the basis of a specified measure of system perform- 
ance (the system evaluation or objective function). 

a * Optimal Design of a Four-Wheeled Vehicle with Direct Mars- 
Earth Communication - C. Pavarini 
Faculty Advisor: Prof. E. J. Smith 

The system model equations (see pages 34-36)' do not include 
volume considerations, for the sake of simplicity. It was found 
that a simple constraint (eqn. 49) could be used to guarantee 
adequate volume for the onboard equipment. 

Based upon a list of equipment priorities determined in light 
of mission goals, approximate relations established for the science 
package are given in equations 1 through 3. 

Analysis of the communications problem revealed that downlink 
(Mars- to- earth) considerations were far dominant. For a pulse-code 
modulated, X-band system assuming worst case link distance, equa- 
tions 4 and 5 relate the important parameters. 

The source of energy for the rover will be radioisotope thermo- 
electric generators (RTGs) . Energy storage will be by nickel-cadmium 
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batteries. Equations 6-13 describe the subsystem. Because 
battery recharge can occur at the end of a Mars day, eqn. 13 
constrains total roving time before a recharge to less than or 
equal to this time period. The square-bracketed term in eqn. 11 
is a probabilistic measure of the percentage of roving time the 
vehicle will be using battery energy. It is based upon exponen- 
tial slope distributions for the Mars surface. 

Temperature control will be achieved through the use of an 
active cooling system with two radiators, on opposite sides of 
the rectangularly-solid equipment package, and an electrical 
heater. Assuming worst case radiation and surface degradation 
effects, eqns. 14 - 30 hold in the steady state. Note that a 
radiator exposed to direct solar radiation will be shut off. 

Based upon analysis of the accuracy of several candidate 
navigation schemes, and relative insensitivity of other subsystem 
operations to navigation error, allocations to the navigation sub- 
system were fixed prior to optimization (eqns. 31 -33). 


The obstacle avoidance subsystem uses a laser rangefinder to 
sense terrain at a distance r a in front of the vehicle, and uses 

this information to calculate approximations of slopes. The term 


2 r AA _ . , 

a ~ refers to the worst case error in slope detection. The 


equation for path length ratio (PLR) was established via Monte 
Carlo simulation. Eqns. 34 - 38 describe the subsystem. 


Computation and data-handling will be accomplished as per 
NASA recommendation, by the use of the Thermoelectric Outer Planet 
Spacecraft (TOPS) computer subsystem, which is described by eqns. 
39 and 40. 


Based on the RPI MRV proposal for a four-wheeled rover, and 
using the fact that the vehicle will be dynamically scaled, vehicle 
structure dimensions and weight can be found by eqns. 41 - 43. 

Finally, equation 44 simply places a limit on the total weight 
of the system to be delivered to Mars. The launch weight used is 
that for a Titan III-C launch vehicle, with allowances left for the 
equipment necessary to achieve a soft landing. 



When the expression that measures system performance is 
written as a function of the design parameters, the optimal design 
problem becomes one of choosing the design parameters to extremize 
that value of the objective function subject to the equality and 
inequality relations of the system model. Since the Mars rover 
has two responsibilities, to rove the surface and perform scientific 
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^estigation, a reasonable objective might be to maximize the 
"uct of distance roved and experimental science time. Of 
rse, there are other possible objective functions. On a Mars- 
Ay basis, in terms of the model parameters, the chosen express- 
i is written: 


T . S . v - T V //c . 

escx sci f rov (45) 

PLR r T ' S . v- T + T +T I 
L sci sex f rov rov rj 

t its of Optimization 

The nonlinear programming problem (NLP) is: 

max f(x) 

F 

F : g^x) ^ 0 i-1,2, . . . ,m ; x^ 0 

x is an n-dimensional vector and the f and g functions are 
aieral nonlinear functions of the elements of x. For solution 
As problem, a program implementing Fiacco and McCormick 1 s Se- 
:ial Unconstrained Minimization Technique was utilized. For 
~LP problem, only local maxima can be computationally located, 

. everal diverse starting points were used to achieve some de- 
of confidence in the solution. In all cases, convergence was 
~he . same point. 

2 

Letting the elements of x be W . , D , P , s*, V-, A, E, , 

sex com 5 com 9 f 9 batt 

and L^., and using the equality relations of the model to 

nate variables, the system design optimization problem reduces 
NLP problem in ten variables (n=10) with thirty-two inequality 
'Taints (m=32). The large number of inequalities results from 
act that many of the eliminated variables must be constrained 
i ive since they represent real physical quantities. 

A first solution of this problem indicated that there were in- 

cient constraints on several system variables due to factors 

• considered in the model. D 2 , s* and V- were too large, while 

com r 

and A were too small. Adding constraints on these variables: 


s* 

4 

20 degrees 


(46) 

v 

g 


1.5 m/sec 


(47.) 

com 


2.0 m 


(48) 

A 

> 

7.35 m 2 


(49) 


\nal solution was: 
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W 


sci 

D 2 
com 


com 

s* 




att 


w. 


L. 

1 


SCI 


66 . 3 kg 
2 m 2 

180.0 watts 
20 degrees 
1.5 m/sec 
7.44 m 2 
9.36 watt-hrs 
3.10 m 
0.0260 m 
0.27 x 10“ 3 /m 


from which all other system parameters can be calculated. The 
constraints of equations 39, 43 - 49 were equalities at the solu- 
tion point. 


b. Optimal Design of a Four-wheeled Vehicle Communicating via 
a Mars Orbiter - F. Stulen 
Faculty Advisor: Prof. E. J. Smith 

Analysis of this problem in the past established that the 
availability of an orbiter would effectively remove all constraints 
on data rate given a 10 kg weight allocation for the communications 
subsystem. Identification of the optimal design in this case then 
requires only that certain modifications are made in the system 
model presented in B.l.a and the NLP problem is resolved for this 
new model. 

Modifications 


The following changes can be made in the system model (new 
equation numbers relate to the unprimed numbers in B.l.a.): 


a. 


The time required to transmit science data becomes 
negligible so that (2) becomes 


T = T 
sci esci 


( 2 ') 


b. Weight allocation to communications is now constant 


W = 10 kg 
com 


(4’) 


c. (5) can be eliminated. 


d. The new optimal design problem reduces to a 9 variable 

problem (as opposed to 10 in B.l.b.) because P can 

com 

be dropped as an independent variable and instead P 

com 

is set to 40 watts. 
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Notation (values of constants are given in parentheses) 


A 

a 


alb(r) 

^alb(s) 


TR 


sun(r) 

^sun(s) 


12 


Dm 


■^batt 

E f 

E fd 

V 

h 


K 

q 

L. 

1 

L 

v 

M 

r 

P 

com 

P 

cp 

PLR 


mv 


P- 


nav 


oa 


prop 

? RTG 

D 

scia 

3 

str 


e 

Si 

Q 


id 


equipment package surface area ,.m 2 . • ■ . ' 

radiometric albedo of Mars, (0.295) 

radiator surface area illuminated by albedo radiation, r? 
equipment package n " " " " 

p 

total radiator surface area, m 

surface area of the top and two sides of equipment package, m 2 
A t plus- area of one radiator, m 2 

radiator surface area illuminated by direct solar radiation, m 2 
equipment package " " " " " y : .» 

ratio of radiator length to height, .(3.33) 

ratio of length of open side of equipment package to its height, (4.0) 
S 1 + a 2 + a i a 2 3 (20.65) 

c ommunicat ions antenna diameter (parabolic dish), m 
battery energy, watt-hrs 

depth of battery discharge divided by efficiency of recharge process, (0 
depth of battery discharge, (0.4o) 

2 ■- 
Mars gravitational constant, (3.75 m/sec ) 

Mars forced convection transfer coefficient, (1.30 watt s/m 2 - °K) 

heat transfer coefficient of insulation, (0.0216 watts/m- C K) 

cooling capacitor, watts/°K • 

insulation thickness, m 

maximum launch weight, (570 kg) 

mass of total vehicle system, kg-mass 

power required by communications subsystem, watts 

" computation and data-handling subsystem, (90 vaccs) 
path length ratio, ratio of actual to straight-line distance roved 
power required by onboard subsystems during roving periods, watts 
power required by navigation subsystem, (6 watts) 

” " obstacle avoidance subsystem, (15 watts) 

power required to propel vehicle, watts 
constant power output of RTGs, watts 
average power required by science package, (26 watts) 
power required by onboard subsystems during recharge node, watts 
" " " thermal control subsystem, watts 

heater power required at night, .watts 
internal heat dissipation during the day, watts 
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Q. 


m 


R 

com 

S 

c 

s* 

S . 

SCI 

T 

t . 

T 

bd 

‘t, 

bn 

T 

esci 

T. 

hi 

T 

intd 


intr 

‘low 


T. 


T 


T 


rd 


rn 


rov 


SCI 


T 

i 

T 
• i 

T 

£ 

•v 

v f 

v t 

w 

c 

W 


com 


cp 


nav 


oa 


SCI 


W 

r 

W 

c 

V? 

I 

W 

c 

W 

w. 


a 

r 

a . 
n 

a 


a . 

si 


internal heat dissipation at night, watts • 

horizontal distance from vehicle to terrain being sensed, ( 30 m) 
communications data rate, bits/sec 

p 

Mars solar constant, (720 watts/ra ) 

maximum slope for allowable vehicle traverse, degrees 

science stops per meter of travel, 

percentage of terrain vehicle will consider impassable 
vehicle track, side-to-side distance between wheels, m 
temperature of. equipment package external skin during day, °K 

" " " " at night, °K 

time required for experimental science during one science stop, sec 
worst case (high) Mars ambient temperature, (275°K) 
required internal temperature during day, ( 310 °K) 

" " " " T}-irrV.+ fPO0°Kl 

worst case (lea;) Mars ambient temperature, (175 °k) 
time required to recharge batteries, hr 
radiator surface temperature during day, °K 

” ** <1 It • 14* O 

night, K 

maximum roving time between battery recharges, hr 

time required to obtain and transmit science data per science stop, sec 
time during Mars day when earth is visible from rover, (10 hrs) 
roving velocity, m/sec 

wheelbase, front-to-back distance between wheels, m 
earth weight of communications subsystem, kg 

" " " computation and data -handling subsystem, (46.8 kg) 

" " navigation subsystem, (15 kg) 

" " ” obstacle avoidance subsystem, (5 kg) - 

" " power subsystem, kg 

" " science package, kg 

" . " " vehicle structure, kg 

" " " thermal control subsystem, kg 

absorbtivity of radiator to solar radiation, (0.5) 

" . " " n infrared , (0.8) 

" " equipment package to solar radiation, ( 0 . 5 ) 

" " infrared " , (0.8) 



AP error in detection of local vertical by navigation subsystem, (0.25 ces) 

e r emissivity of radiator, (0.8) 

G s emissivity of equipment package skin, (0.8) 

l- 1 ^ coefficient of kinetic friction of vehicle, (0.10) 

<T Stephan- Boltzman constant, (5.67 x 10 - ^ vatts/m^-°K*) 
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Optimization 

The changes noted above will obviously lead to an optimal 
design quite different than that identified in B.l.a. This new 
design will indicate the gains in system performance to be de- 
rived from servicing the MRV with a Mars-orbiter which has the 
capability of relaying data from the vehicle to earth. 

Identification of the optimal design has not been completed, 

; hut is in the programming stage. Results are expected by March 15, 
1973. 

B.l.c. Optimal Design of a Six-Wheeled Vehicle with Direct Mars- 

Earth Communication - M. O'Callaghan 

Faculty Advisor: - Prof. E. J. Smith 

Modifications 

Use of a 6-wheeled vehicle (e.g. Bendix or McDonnell 
Astronautics prototypes) requires considerable modification to 
the original system model. Major differences appear in the thermal 
control, obstacle avoidance, (and obviously) vehicle structure sub- 
systems. 

The six-wheeled vehicle consists of three segments, each with 
two wheels, coupled together. The thermal control problem is much 
different from the 4-wheeled case where there is only one equipment 
package to control. For the 6-wheeled configuration, there are two 
separate volumes requiring temperature control, each with different 
internal heat dissipation. Thus, the magnitude of the thermal con- 
trol problem doubles, and in addition the two new problems are 
coupled. 

The mobility of the 6-wheeled vehicle is quite different from 
the 4-wheeled version. Thus, (38) no longer holds, and a new ex- 
pression is necessary. 

Finally, the equations of the vehicle structure itself must 
change. The most important modification occurs in (43), which re- 
presents the payload fraction consideration. While for a 4-wheeled 
vehicle it is expected that the ratio of equipment weight to 
structure weight can be as high as 2:1, the 6-wheeled case should 
decrease this to 1.75:1. 

All of the required modifications to the system model have 
been completed. 

Optimization 

As in B.l.b., the identification of the optimal design has 
proceeded to the programming stage. However, since this problem 
is more complex than that of B.I.4., more problems can be ex- 
pected in obtaining the NLP solution, and so the projection for 
completion in this category is April 15, 1973. 
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B.2 Time Simulation of Vehicle Performance - R. Gallerie and T. Mans 

Faculty Advisor: Prof. E. J. Smith 

This is a new task area developed for the purpose of investi- 
gating the time performance of the vehicle. There are two motiva- 
tions: 

1. The optimal designs were established from static 
models which accounted for dynamic processes by 
averaging and/or worst-case approaches. It is 
important to check the validity of these approaches 
by finding out exactly what will occur onboard the 
vehicle during its operational life. 

2. Dr. Garrett Paine of JPL has indicated interest in 
how the vehicle would perform under less autonomous 
operation, i.e. if earth commands were more frequent 
and specific. A time simulation will allow just 
that kind of analysis. 

The goal in this area is to develop a working simulation of 
the vehicle by 1 June 1973. Major considerations will be flexibility 
of the simulation package, complete description of the vehicle time- 
state, and capability to simulate a comprehensive set of stochastic 
inputs and mission profiles. 

B.3. Mars Roving Vehicle Computer - C. Boroughs and F. Samuel 

Faculty Advisor: Prof. E. J. Smith 

The objective of this task is to determine .the onboard computer 
requirements for the Mars Roving Vehicle and to translate these re- 
quirements into a practical device that can be built and tested with 
the interfacing vehicle subsystems. 

The onboard computer must be a very sophisticated device. Be- 
cause of the communication delay time between earth and Mars and the 
nature of the mission profile, the computer must have the ability to 
make decisions concerning path selection, obstacle avoidance and 
navigation. It must also monitor, sequence and record the results 
of the MRV's scientific investigations as well as monitoring and 
controlling its own internal housekeeping functions which include 
power and temperature control and earth communication. 

A study of the subsystem functions with respect to computer 
design has been initiated. Mathematical models of subsystems are 
being revised to be adaptable to the computer design. The computer 
will be modeled on a queing theory basis. The machine queing will 
be dictated by the subsystem functions. Then various computer 
architecture can be compared using the queing theory pattern de- 
veloped. 

In addition to definition and evaluation of parameters for the 
subsystems, preliminary design work was done in support of the 
vehicle group, Task A. The control program and hardware interface 
for the GEPAC-30 have been designed. This will be used to control 
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the 0.4 scale model under development. Effort has been directed 
at the design of an interface between the proposed hand held 
radio control unit and the existing hardware on the GEPAC com- 
puter. 

In the next period, the complete subsystem parameter evalua- 
tion with the queing model will be completed. The design, writing 
and debugging of this simulation program will be effected. 

Work on the 0.4 scale model and its associated electronics 
will be completed. It is anticipated that a certain amount of in- 
structional effort will be required to acquaint the vehicle group 
with the operation of the GEPAC-30 computer. It is anticipated 
that computer programming will be required to support this effort. 

Task C. Navigation, Terrain Modeling and Path Selection 

The mission plan to undertake a systematic exploration of Mars 
requires that the roving vehicle can be instructed to proceed under remote 
control from its landing site to a succession of desired locations. This 
objective requires that the vehicle possess the capabilities : of sensing 

and interpreting the terrain to provide the information required by a path 
selection system, of selecting paths with due regard to safety and other 
considerations, and of knowing its location and that of its destination. 
Emphasis has been directed to the subtasks of: range measurement using 

lasers, obstacle detection systems, terrain modeling, and path selection 
systems evaluation. Details of the progress made are provided below. 

a * Design of a Laser Rangefinder for Terrain Modeling - D. Palumbo 
Faculty Advisor: Prof. C. N. Shen 

The objective of this task is to investigate alternative 
methods by which range may be determined using a laser system. 
Included in the study were: optical focussing, timed pulse and 

phase difference methods. A hybrid scheme combining the timed 
pulse and phase difference concepts has been defined. 

Optical Focussing 

This method was investigated primarily for use with the 
"Stereo Angle" scheme being investigated under Task C.l.c. and 
involves the following basic features: 

1. Laser light is used to illuminate a point. 

2. A detector scans the field, "looking" for the 
point of laser light. 

3. As the detector approaches the point, the in- 
tensity the detector "sees" will begin to rise 
until a maximum is reached and passed at detector 
angle 

If data was stored during the scan concerning intensity vs 0 , 
the value £)' at maximum intensity can be determined through com- 
puter methods. Once ©' is known, geometry yields 

tan Q‘= R/B 

where R is the range and B the laser- detector separation. 


u 
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The accuracy of this method depends on how well 0 can 
be resolved. Calculations made under Task C.l.c. indicate that 
accuracies on the order of five arc seconds are necessary to de- 
termine R within five centimeters ^ Although this may at first 
seem formidable, encoders are available at this accuracy. Gurley- 
Teledyne of Troy, New York, produces such a mechanism. It consists 
of an encoder wheel divided with 16,384 slots per revolution. The 
slots are sensed by a LED-photodiode combination. Using electronic 
extrapolations, accuracies of 4-5 arcseconds are available. The 
encoder is cylindrical, 6 inches long with a 3 inch diameter; 
electronics are rectangular with dimensions, 3x4. 5x2 inches. The 
entire system weighs 5.5 lbs. and can be bought for 4-5 thousand 
dollars. Output from the encoder is in the form of a digital up- 
down counter. 

Although the system seems feasible at this point, many factors 
weigh against it. The one dominating flaw is a resulting slow scan 
rate * A fast scan rate would eliminate error in data acquisition due 
to vehicle rock and roll. The necessity of the system to scan for a 
point, plus added encoder inertia would limit performance beyond 
fast scan assumptions. 

Timed Pulse 

In the Timed Pulse method, the laser is pulsed and split into 
two beams. One represents the leaving beam, the other the returning 
beam. The leaving beam is used to start a counter, which records 
elapsed time of flight, dt, when stopped by the returning beam. This 
time is related to the range through the speed of light. 

In a previous investigation of the Timed Pulse method at Rens- 
selaer , it was found that in order to achieve a 5 cm. accuracy 

a time resolution of .25 nanoseconds must be available for dt. Ob- 
taining this resolution is a formidable task. For example, Reference 
6, in 1969 investigated obstacle detection for the lunar rover, and 
reported accuracies of 30 cms. In 1971, W. Kuriger, Reference 
states that under laboratory conditions, accuracies of 10 cms. were 
available. Further investigation showed that a counter is com- 
mercially available through Nanofast corporation. The counter has 
.25 nanosecond resolution, but comes in large modular form, is heavy 
and lists from 4-5 thousand dollars. 

Phase Difference 

The basics are as follows. If the laser beam is modulated at a 
frequency f, the resulting wavelength, L, will be given by 

L = c/f 

where c is the speed of light. The rangefinder is to have a depth 
of field R in meters such that, 

10 £ R ^ 50. 


The maximum distance the wave will travel will then be 100m. So 
L must be greater than 100 meters to insure a unique relation be- 
tween the phase and the range. At L = 100 m the frequency of 
modulation will be given by 

f - c/100 3Mhz 

A pulse generator will be used to turn the laser on and off. 
When "on", a pulse modulator will superimpose the megahertz signal 
on the pulse. The modulating amplitude should not be large enough 
to turn the laser on or off. If the threshold voltage of the laser 
is given by Vt, then the modulating amplitude should be approxi- 
mately equal Vt at most, so that the pulse amplitude must be twice 
Vt. This allows maximum modulation without disturbing laser opera- 
tion. The laser pulse will contain several cycles of modulation. 

The phase difference between the reference, Vr, and the detector 
voltage, Vd can be determined through a heterodyning technique, ex- 
plained below, (see Figure 11). 

When two signals of different frequencies (wl, w2) are com- 
bined, side band frequencies are generated at the sum and difference 
frequences. If 


• Vd = sin (w^t + 9 ) 

Vm = sin (w^t) 

where Vm is the modulating voltage, then, 

V = sin (w^t + 0) - 1/2 cos (Wj+W£+ 9) 

If Vd and Vm are cosine functions then 

V = cos (w^t+9) +1/2 cos (w^t+W 2 t +&) 

1/2 cos(w^t+W2t+0) 

These relations can be proved through Euler's equations 

exp(iy) = cos(y) + isin(y) 

exp(iy) + exp(-iy) = cos(y) + isin(^) + cos(y) 
- isin(y) 

cosy = (exp(iy) + exp(-iy))/2 
cos(x)cos(y) = l/4(exp(i(x+y))+ exp(i(y-x)) 
exp (- i (y+x ) ) +exp (- i (y-x ) ) ) 

= l/2cos(y+x)+l/2cos(y-x) 

If we let y = w^t+ d 

x = w 2 t 

the relation is proved. 

Using proper modulation frequency and band pass filtering, Vd may 
be reduced from 3Mhz to 3Khz and compared with Vr at 3Khz. 


ULSE PULSE 



Figure 11. Basic Phase Difference Technique 




If for example 0 = "tt/2, this corresponds to one quarter of 
•the wave*s period at 3Mhz or 83.3nsec. A phase delay of IT/ 2 
corresponds to a range of 25 meters. The resolution needed for 
5cm. accuracy is 

. 05/25(83. 3nsec) = .16nsec. 

At 3Khz. , after heterodyning, the period will be one thousand 
times greater, 83.3 microsec. and .16 microsecond resolution will 
be needed to resolve five centimers. 

Conceptual Design Based on Hybrid Method 

The Hybrid Method combines the Timed Pulse and Phase Differ- 
ence methods to achieve accuracies within the requirements neces- 
sary. This type of design is used in surveying equipment with 
accuracies reported at 1 cm. 

If a range calculation has just been completed and a new 
range is to be determined at a different point, the pulse generator 
is signaled. The pulse turns on the laser allowing modulation by 
the frequency generator. The returning wave is reduced by mixing 
with appropriate frequencies produced by the frequency generator 
until a 3Khz wave results. Zero detectors are used for controlling 
the counter by sensing the reference and reflected signals. The 
time displayed on the counter will represent the elapsed time of 
flight. Since many zero crossings can be obtained per pulse, the 
range may be averaged to further improve resolution. At the end of 
the calculation the rangefinder is ready for another pulse, (see 
Figure 12 } for details. 

Another design consideration is if the reflected beam will 
have enough power. Although the device will be operated primarily 
in the laboratory, some insight might be gained from looking at 
W. Kuriger’s calculations for a Martian Rangefinder, Ref. 6 

Pr = received power 

Pt = transmitted power = 15 watts 

E = optics efficiency = .4 

Z = Martian surface reflectivity = .3 

Ar = area of receiver = .49xl0”2 m 2 

then, 

Pr = PtEZAr/2piR 2 

For a range of R = 30 meters, Pr equals .137 microwatts. At 5 
meters, Pr equals 5 microwatts. With this in mind the following 
components were chosen for the prototype. 

GaAs single diode laser 

RCS TA 7609 @ 15 watts 

photodiode - UDT Pin-5 2 

Active area - .037 cm. 

NEP - 6.5xlO~^w/Khz 

- 9 

Risetime - 5x10 sec. 
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Figure 12. Final Design 




Timetable for Future- Investigation 

FEBRUARY MARCH - APRIL MAY 

laser operation —————— — 

analysis of 
. terrain reflectance 

incorporation of 

heterodyning 

conclusion — 

C.l.b. Terrain Sensor Instrumentation- Laser Scan Methods - G. Herb 
Faculty Advisor: Prof. C. N. Shen 

It is assumed that a laser range finder type of sensor will 
be used to supply the terrain data used by the Mars roving vehicle 
during its autonomous navigation. If this data is to be collected 
as the vehicle is moving, the rocking and rolling of the vehicle 
as it travels the uneven terrain can be a large source of error. 

The objective of this task is to reduce this error. Accordingly, 
the laser scan methods considered here are of a block scanning 
nature. That is, the entire field of view of the rover is scanned 
out very rapidly so that the rock and roll error between points is 
negligible and all points in the field may be considered to have 
been obtained simultaneously. 

The criterion for a system worked under were: 

1. No or few moving parts 

2. A minimum of 10 data points per horizontal scan 

3. An angle incremented scan (i.e. point spacings 
vary with distance from the vehicle) 

4. Scan requirements: azimuthal angle 40° - 60° 

range angle - 40° 

The last criterion results from the projected dimensions of 
the vehicle and the range of the scan (about 3 to 30 meters). The 
RPI vehicle is currently designed to be 10 feet wide, or about 3.1 
meters. A 60° scan angle (i.e. + 30° from the direction of travel) 
is required to cover this width with a scan at 3 meters from the 
vehicle whereas at 5 meters a 40° scan angle is required. 

The first method investigated was a system using oscillating 
mirrors to do the sweeping. With this system it is very easy to 
synchronize the transmitter and the detector since both use the 
same optics. The transmitted pulse is fired, and the return pulse 
comes back through the optics a few nano-seconds later. The 
mirrors have essentially not moved during this time and thus the 


detector is "looking" at the same spot at which the pulse was fired. 

The problem with this system is that it uses oscillating 
mirrors from which the angle data must be taken. At the scan fre- 
quencies of interest (> 10 KHz) angle decoders cannot be used on 
the oscillating mirrors because of the decoder disc's non-negligible 
inertia at those frequencies. Spinning poly-faced mirrors would 
solve the inertia problem, but synchronizing the two spinning mirrors 
with a high degree of accuracy and keeping them synchronized would be 
an extremely difficult task. It is therefore felt that this system 
. would not be satisfactory. ■ ^ 

Electro- Optics 

In order to decrease or eliminate this necessity of detecting 
angle information from a mechanically oscillating mirror, it seemed 
that a reasonable approach would be to investigate electro-optic 
and other electronic scanners. 

There are several types of electro-optic devices. The three 
considered were: (a) acusto-optic beam benders, (b) piezoelectric 
deflectors and (c) digital light deflectors. 

Acusto-optic beam benders use the property of some crystals 
which causes their refractiveindex, n, to have a gradient propor- 
tional to a pressure gradient set up in the crystal. Therefore, 
since light is bent toward an increasing gradient in the refractive 
index, by changing the index gradient of the crystal the deflection 
of the beam is changed. In crystals which possess the acusto-optic 
property, the index gradient may be changed simply by changing the 
pressure gradient in the material. The major problem with this 
type of deflection, however, is that the proportionality constant 
between the pressure gradient and the n-gradient is extremely small, 
Ref. -7. Typical units provide deflection angles on the order of 
tenths of degrees. These angles are much too small for our purposes. 

Piezoelectric deflectors use the well known piezoelectric effect 
of certain crystals which causes them to deform in certain ways when 
an electric potential is applied across them. The most common de- 
formation used are the longitudinal and the shear deformation. A 
reflective surface is attached to one face of the crystal, and it is 
the motion of this surface as the crystal is deformed which is used 
to cause the deflection. 

The advantages of this type of deflection is that the trans- 
mitted light pulse and return pulse may use the same optics. It 
thus has all the advantages of the vibrating mirror scan described 
earlier without the problems associated with the mirrors being . 
mechanically scanned and synchronized. But once again max imum angle 
is limited. Commercially available units have maximum scans on the 
order of 6°, Ref. 8. Also the angle data must be obtained indirectly 
from the analog control voltages. 
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Digital deflectors are sensitive to the incident light’s 
direction of polarization. The incident beam will be laterally 
displaced or angularly deflected to one of two possible positions 
depending on its direction of polarization. These deflectors 
consist of a polarization modulator, such as a Pockel or Kerr cell, 
which either rotates the direction of polarization 90° or allows it 
to pass unrotated, followed by a binary birefr ingen t discriminator, 
i.e* a crystal with two indices of refraction depending on the di- 
rection of polarization of the incident light. Much work has been 
done on these deflection devices, notably by Kulcke, et al, Ref. 9 
and T. J. Nelson, Ref. 10. The number of positions available at 
the output of a deflection system can be increased by cascading 
individual binary units. Properly cascading N of these units re- 
sults in 2 n positions in the final output. Going from a one 
dimensional to a two dimensional scan is easily accomplished by 
adding an M unit scanner rotated by 90° with respect to the first 
N units. This results in an array of 2^ by 2^ possible deflection 
positions. 

The binary units accomplish a lateral displacement as shown 
in Figure 13a* An incident beam of collimated light will be deflec- 
ted by an angle O if its polarization direction is in the plane of 
the crystal’s optical axis. If the direction of polarization is 
normal to the plane. of the axis, the beam is not deflected at all. 
When the deflected beam hits the edge of the crystal, it is again 
deflected so that it emerges from the unit parallel to the incident 
beam. Therefore a two dimensional digital scanner will provide a 
2^ by 2^ matrix of parallel beam positions. A method of converting 
these parallel beam scans into angular scans will be covered later. 

The angle 0 is a function of the crystal material and the 
angle the optic axis has with respect to the incident beam. The 
figure also gives maximum values of d for three crystal materials 
and the value of ^ needed to obtain that maximum. 

Note that in each cell the control voltage is either present, 
or absent, corresponding to either a rotated polarization in that 
cell or a non-rotated polarization. Therefore, each position 
available at the output can be associated with a unique binary code 
word determined by which cells are energized and which are not for 
that point. Thus each point in the sweep is known exactly by the 
particular code word being issued by the control logic at that in- 
stant. 

There are several problems with these binary units. Two are 
mentioned here for completeness. 

To rotate the polarization of light by 90° in a Pockles cell, 
the half wave voltage must be applied to the cell in such a way as 
to set up an electric field in the cell parallel to the direction 
of light propagation. To accomplish this, electrodes must be 
placed in the optical path. This then creates the problem of 
finding the optimum compromise between the two competing requirements 


of low resistance electrodes and maximum transmittance of light 
by the electrodes* Kulcke found highly acceptable results could 
be obtained using a variety of electrode-adhesive combinations, 

Ref* 9. It was also found that the modulator electrode trans- 
mittance was the limiting factor and hence the total system 
transmittance could be approximated as 

Tsys = t 2n 

where t is the transmittance of an electrode- adhesive combination, 
and n is the number of cells in the system (2n appears because 
there are two electrodes per modulator). 

Another problem of discriminator cells is background light, 
that is light appearing at a point or points other than the one 
called for. The major cause of this is deviation of the modula- 
tors' control voltage from the required half wave voltage. As the 
control voltage becomes greater than or less than the half wave 
voltage, the polarization of the light is rotated more or less than 
90°. It will then have a component in the unrotated direction. The 
two components will be split by the discriminator crystal and the 
undesired component will appear as ghosting at the un-called-for 
position. Values of this background light have been graphed, 

Rsf. 9, and are shown in Figure 13b. Note that the background 
light can be kept beJLow the incident light by more than 40 db if 
the voltage deviation is kept below 10%, a tolerance which is 
easily obtained with a good regulator circuit. 


The advantage of this system is that the positions and there- 
fore the scan angles transmitted are known exactly. A position is 
called for by the scan logic and that scan angle (both horizontal 
and vertical angles) is produced. Therefore no indirect measure- 
ment is needed, and the associated error is thus eliminated. The 
disadvantages are its relative bulk, and the large voltages re- 
quired (the half wave voltage of a Pockels cell is on the order 
of a kilovolt). Also the return beam cannot use the same optics. 

An alternative to the digital deflector optics is a matrix 
of N by M individual laser diodes. This has all the advantages 
of the binary deflectors, but lacks the problem of high voltages 
and is somewhat less bulky. It does, however, require N times M 
diodes whose alignment is critical. 


Converting Parallel Positions to Angies 


The last two methods both provide parallel light outputs. 
This must be converted to scan angles if the method is to be use- 
ful in a scanning system. 


This can be accomplished by placing the origin of the parallel 
positions on the focal plane of a small apature lens. This is 
shown in Figure 17 for the case of a matrix of laser diodes, i.e. 
uncollimated light. The case of binary deflection optics with 
collimated light is similar except that instead of collimating, 
the lens simply deflects the beam almost as if it were an ideal 
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ray". 
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Scanning Systems 

One possible system, the dual vibrating mirror scan, was 
mentioned much earlier in this task summary. Here, two more 
systems will be introduced which attempt to alleviate the problems 
of the dual mirror approach. 

Kuriger, Ref. 11, has worked on and proposed a hybrid scan 
system using a vibrating mirror for a vertical scan and a line of 
laser diodes to accomplish the horizontal scan. His range finder 
system is not capable of producing the necessary accuracy, but his 
scanning system does have some advantages. Kuriger 's system is 
presented with some minor modifications in Figure 15. The "pulse 
transmitted detector" has been added to conform with the information 
needed in range finding schemes currently under study. In this 
system both the transmitter and the detector use the same vertical 
scan mirror. The diodes provide the horizontal scan, while the 
receiver looks at the entire horizontal scan area via the cylindri- 
cal lens. The use of the same scan mirror assures synchronization 
between transmitter and receiver, and the narrowing of the detector's 
field of view to only one horizontal row at a time increases the 
signal to noise ratio. Further advantages of this system are its 
compactness and small number of components. But although the hori- 
zontal transmitted angle is known exactly, the vertical angle deter- 
mination relies on the same type of mirror position detection as 
does the dual mirror system and suffers the same problems. 

An alternative system which seems to avoid all the major scan 
problems is now presented. Figure 16 shows the proposed system 
using digital deflection optics, but a matrix of laser diodes could 
be substituted just as well. 

The detection system used is similar to Kuriger 's except that 
the detector section is mechanically a separate entity from the 
transmitter. When the mirror drive starts its vertical sweep the 
synch (start) circuit detects this and initiates the transmitter's 
scan. The transmitter's vertical scan is approximately synchronized 
to the mirror, but the detector is given a wide enough field of view 
to compensate for a slight mis-match. Any mis-match which does 
occur does not continue to accumulate since the transmitter and de- 
tector are re- synchronized at the beginning of each vertical scan. 

The advantages of this system are many- fold. First, the 
actual laser scanning is completely solid state. Also the trans- 
mitted angles are known exactly with each pair of horizontal and 
vertical scans having its own binary code lable. Finally, although a 
hybrid detector scheme is used, data point accuracy is extremely im- 
proved over Kuriger 's scheme because the mirror position is not used 
to give point location but only to narrow the detector's field of 
view to the approximate azimuthal scan line of the transmitter. This 
system then seems to have the potential for extremely accurate 
transmitted angle information while maintaining a reasonable signal 
to noise ratio at the detector by narrowing its field of view. 

The maximum scan angle obtainable with this system is determined 


by the angle generating lens system used. For the single lens 
Set up shown in Figure 14, maximum scans of 55° can be obtained 
with a good lens. If a multiple lens system' is used, this 
Viewing angle" can be greatly extended from this maximum. For 
example, using standard photographic lenses as a guide, a 28 mm 
focal length lens (6 elements) has a viewing angle of 74°, a 20 mm 
lens (11 elements) has a viewing angle of 94°. If distortion of 
field is allowed (an effect which can be compensated for by appro- 
priate placement of beam sources on the matrix plane) fisheye 
lenses can also be used. Fisheye lenses have viewing angles from 
.180° (10mm focal length) to 220° (6mm focal length). Reference 12. 

Since all lenses are circular, the angle maximum is completely 
symmetric with respect to the optic axis and therefore holds* for 
both azimuthal and range angles* 

Scan rate is limited only by position switching time. For 
the case of a matrix of laser diodes, minimum switching time is 
caused by the diodes 1 minimum pulse width and the maximum usable 
frequency of the driver circuits. This is on the order of magni- 
tude of 100 KHz or 10 micro seconds per position. Thus the limiting 
factor for sweep rate will not be the scan unit, but rather the mini- 
mum on-time of the laser needed by the range finder system. This 
will be much larger than 10 micro seconds. 

For a matrix of diodes system the weight of this scan mechanism 
will most likely be under 10 pounds. 

One problem with all the angular scan methods is the incident 
spot distortion at the terrain. Since a collimated beam is trans- 
mitted at an angle, the actual shape and size of the spot will be 
distorted when it strikes all surfaces other than those to which it 
is normal (see Figure 17). This is most severe at the maximum 
ranges. As can be seen from' the figure 

w t _ W for horizontal terrain (i.e. in 
sin 0 the plane of the rover) 

Therefore at 30 meters (0 = 5.7°) a 2.5 cm beam diameter ( = 1 inch) 
would produce a spot length of 25 cm. This effect worsens for nega- 
tive in path slopes and improves up to a point for positive in-path 
slopes. 

This spot distortion phenomenon will be analyzed in greater 
detail later with respect to its resulting phase distortion or 
blurring of the modulation frequency of the laser. 

Shown below is the projected schedule of work to be undertaken. 
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Final Report - 

C.l.c. Stochastic Estimates of Gradient from Laser Measurements - P. Burger 
• Faculty Advisor : Prof. C. N. Shen " 

The objective of this task is to determine the errors in 
estimating gradients from laser range measuring systems as a result 
of measurement errors for the purpose of path selection decision 
making. Autonomous terrain sensing systems of high reliability will 
be required by an unmanned mission because of the round trip com- 
. munication time delay. At this time, the investigation is focussed 
on determining terrain features in the 3 to 30 meter range. 

The central concept of the system is the range finder, which 
locates a point on the terrain according to a range measurement R, 
azimuth angle @ , and elevation angle ft . Two points along the 
path of the vehicle determine an in-path slope, while those across 
the path compute a cross path slope. However, inaccuracy in measur- 
ment can introduce very large errors in the computed slopes and 
heights, which are the main factors in path selections. There are 
some threshold values for these factors above which a change of path 
is required. Ref. 13. 

In brief*, the results which have been obtained to date relate 
the standard deviation of the gradient which an operating system 
would encounter as a function of the standard deviations in the 
range, deviation and azimuthal measurements of a scanning system. 
Accordingly, the consequences of accepting hardware laser systems 
with specified accuracy limits can be assessed in terms of the un- 
certainty in the calculated gradient and ultimately in the probabi- 
lity that an improper path selection will be made. These results 
also provide a basis for specifying optimum data spacing. 

Method of Approach 

From the measurement data it is desirable to obtain the maximum 
slopes and the elevations of the terrain in front of the vehicle. 

a. Transformation of Coordinate Systems 

The quantities R, 0, and jS are measured with 
respect to the coordinate system h", a M , and b", fixed 
to the vehicle, Figure 18. With laser height fixed at 
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3 meters : 


h" = 3 - Rsin (3 
a" = Rcos /S sin & 
b' = Rcos/3 cos O 

The body-bound axis rolls with the angle and pitches 
through an angle £ about a reference frame h, a, and b 
formed by the local vertical and an axis in- a plane con- 
taining the heading and the local vertical. The coordinate 
transformation, Ref. 14, is: 
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= C(*rf) B(p 


h" 

a" 

b" 


( 2 ) 


where C(^) and B(|) are the roll and pitch transforma- 
tion matrices respectively. 


b. Determination of Slopes and Gradient 

A number of measurement points within a certain 
small area of surface can be used to determine a plane 
in space, which may be written as: 

h = ax 1 + bx 2 + x 3 (3) 

where is the cross path slope and X 2 is the in-path 
slope, which are the two constant parameters to be 
determined, Ref. 15 and 16. The corresponding gradient 
of the plane, which is defined as the gradient of the 
terrain is 

Sg = . (X 1 + X 2^ (4) 

If the slope is less than a predetermined criterion, it 
is considered to be safe for the vehicle to travel ahead. 

In order to locate the plane by a number of measure- 
ment points, one may rewrite Eq. (3) as 


h. 

x 


a i X i 


+ + x 3 


(5) 


where h^ a^ and b^^ are found from R^ © i , and in 

Eqs. (T) and (2) for each ith point. Theoretically, three 
points determine a plane (i=l,2,3). For greater accuracy 
however, more than three measurements (i=l,2,3, . . ,n) 
(probably n=4 or 6) are needed to determine the slopes. A 
complete picture of the terrain in front of the vehicle can 
be constructed by modeling numerous adjacent planes, each 
covering a small area of surface. 


Analytic Solutions for Parameter Estimation 

In this section are listed the solutions of the least square 
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estimate, the covariance matrices when n ? 3. 

a. Least Square Estimate of Slopes 

If aj = aj and are assumed to be true in 

E<1« (5), a least square error estimate can be performed 
which minimizes 

£ „ _ 2 

it! (h i - h i> 

where hi is the actual measured height and h^ is the 
corresponding height in the modeled plane. 

In matrix -vector notation, Eq. (5) is: 
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n = 4 or 6 


The least square estimate of the parameter x becomes. 
Ref. 17: 

A -T- -1-T- 
x = (A A) A h 


b. Perturbation of the Variables 

The least square estimate in the previous section 
assumes the matrix A to be completely deterministic. 
Ref. 18. In reality, however, there is error involved 
in the determination of a, b, and h, due to the errors 
in the measured values of fi, f , R , f3 , and 0. 

If the symbol S denotes a perturbation, thenbh, 
Sa, and Sb in terms of SR, £/i, and be are: 


c. 


ih 

Sa 

6b 

r — 


D(h", a", b", 



C(0)B(tj)G(R, (3,0 


Covariance Matrix of the Variables 


The covariance matrix of the variables can be de- 
fined as 


M ■ E { 
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where E denotes expected value. 
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If S$y and 56 are not correlated, 

then with the aid of Eq. (9) 
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(11) 
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From Eqs. (10) and (11) the standard deviations of 
h, a, sind b, can be computed in terms of those of 5^ ., 

§ > R > £ > and 0 for each point. These are known 

.quantities that depend upon the accuracy of the measur- 
ing devices. 

d. Covariance Matrix of the Slopes 

Eq. (5) can be written as 

h = Ax (12a) 

If h = h + Sh, A = A + £A, and x = x + Sx, (12b) 

then with the aid of Eq. (7) the covariance matrix of the 
slopes is found as: 


E [S x Sx T J = F ^ E ^h£h T j 

where fe 

U * 3 


E ^AxSh T ] - E[Sh(gSx) T ] + E [£ Ax(SAx) T ]"^ F T ' (15) 



Z A* - 


t Sb,x/ 


(16) 


This equation utilizes the least square estimate of 
the parameter x from Eq. (8) and the elements of the co- 
variance matrix of the variables for each point, which are' 
determined from Eq. (11). 

e. Variance of the Gradient 

If the symbol <5^ denotes the standard deviation of 
the gradient, Sg, then from Eq. (4) 


dSg = (Xj^ + x 2 ) 2 x^dx^ 


-r 



The value of 0^. gives a rough estimate of the 
accuracy in the estimation of the gradient. If the esti- 
mated gradient, Sg, is 20° with <53j- = 2°, then there is a 
68% probability that the actual slope is between 18° and 
22° and a 95% probability, Ref. 19, that it is between 


( 17 ) 


16 and 24 . In general, there is a 68 % probability 
that the actual gradient lies between the values of 
Sg + and Sg - and a 95% probability that the 

actual gradient lies between the values of SG + 2 
and Sg 2^sj. This is important when the estimatidn 
is close to the maximum permissible value. 

Numerical Results / 

It is assumed the rover will use a split beam where A/S is 
the difference in the elevation angle between the two beams and 
A 0 is the difference in azimuth angle between any two laser 
pulses as shown in Figure 19. 

The values of /3 and 0 for each data point, along with the 
magnitude of the cross-path and in- path slopes of the terrain, 
determine the data point spacing and the distance of the data 
points from the vehicle. After a number of data points are 
measured, they are transformed to the non- rotating reference 
system by Eqs. (la to 2c). Then a least square estimate of the 
slopes is performed which is given by Eq. ( 8 ). This least square 
estimate depends directly on the values of h, a, and b for each 
measurement point. Utilizing the value of the vector x f round in 
Eq. ( 8 ), the values of 0 , % , R, (3 , and B for each data point, 
and the standard deviations of these quantities, the magnitude of 
the gradient and its variance can be determined. 

The standard deviation in gradient, 0*3 is calculated by 
utilizing Eqs. (11), (14), (15), and (17). From Eq. (15) and 
(17) it can be noted that the variance of the gradient depends 
in part upon the cross path and in path slopes of the terrain. 

Eqs. (11) - (15) depend ultimately upon the measured values of 
, £ , R, {3 , and 0 for each point, and the standard deviations 
in these quantities. That is to say, the variance of the gradient 
depends additionally’ upon the roll and pitch angle of the vehicle, 
the data point spacing, the distance of the data points from the 
vehicle, and the accuracy of the measuring devices. How much each 
of these factors affects OV^ is to be investigated. 

For simplicity, the roll and pitch of the vehicle are first 
set equal to 0 , and the reference gradient of the measured plane 
is assumed to be 0° wherever it may be located. By setting 
0 & = = 0 . 1 ° and <fy= <3| = 1 ° it is determined that for any 

reasonable data point spacing at 20-30 meters from the vehicle, 

was on the order of 30°-60°. On the other, hand, if <^= Gf = o* 
Cs<^ decreased to only 2°-3° at that distance. Since we can expect 
that <5jjr and a§ are really in the neighborhood of 1° (due to the con- 
stant rock and roll of the vehicle as it traverses the surface), 
the rover will have to be provided with a "rapid scan" laser. If 
the scan rate is on the order of milliseconds, which is perfectly 
feasible with electronic scanning, then each set of 4 adjacent 
data points is measured practically instantaneously, since the 
rover motion is on the order of seconds. Each of the 4 data points 
will have the same value for <f and § and therefore all 4 points 
will retain the same relative position to each other when they are 


59 


transformed from the vehicle coordinate system to the fixed 
system. Consequently, the rover can model the planes in the 
vehicle coordinate system and then transform the planes to the 
non-rotating frame. Since the maximum positive or negative 
slope that the vehicle can navigate is + 25°, relative slopes 
might be as high as + 50° and still be navigable, as shown in 
Figure 20. Therefore, a range of slope changes from +50° to -50° 
have been considered. 

Also due to rapid scan, the error introduced bv setting 

= 1° will be consistently between 1° and 2° because it 
only involves the error in transforming the already modeled plane 
from the vehicle system of reference to the fixed system. Conse- 
quently, the effect' of and is not considered in this analysis. 

On the basis of a flat terrain, with = S/a = 1' and^R, 3 5 cm, 
and vary A/iand b& (and consequently the data point spacing ) for 
each set of 4 points, the graphs shown in Figure 21 are obtained. 

Each solid line represents constant values for A6 and A (5 . Ab and 
A a are the data point spacings along the in- path and cross path 
directions respectively. It is noted that for any constant value 
of LB and &£> , the data point spacing decreases very rapidly as the 
scan approaches to the vehicle and consequently the error in gradient 
rises very rapidly. 

In choosing an "optimum" data point spacing it must be kept 
in mind that by increasing the spacing the error in gradient is 
decreased. However, this renders the data less meaningful, as more 
terrain has been overlooked. At the 30 meter range, a general 
picture of the terrain with a spacing of 2-3 meters is sufficient. 

At close range, (4-7 meters) the points should be at least as close 
as 0.66 meters because this is the width of the widest navigable 
crevice. In Figure 21, the dotted line represents ^ vs distance 
from the vehicle for such a scheme, where the data point spacing 
varies from about 0.5 meters at a distance of 4 meters from the 
vehicle to about 3 meters at a distance of 30 meters from the vehicle. 

Each of the plots in Figures 22 and 23 utilize the same optimum 
scheme as the dotted line in Figure 21. Once again flat terrain is 
assumed for all the data points. In Figure 22, is kept at 5 cm 
and gander© are set to 0, .017 (1‘), and .1 (6'). The error in 
gradient is markedly reduced by decreasing <r?a and cr e from 6' to 1' 
but any further increases in accuracy will result in bulkier and 
heavier equipment which is not very beneficial. By comparing the 
plots for = °© = 1' and = (T& = 0° it is seen that a standard, 
deviation in <T a and cf e of 1' has very little effect on the error in 
gradient. 

In Figure 23, is varied from 1 cm to 10 cm while and G& 
are kept constant at 1' (.017). The error in gradient increases 
rapidly at close range for each value of Or . 

Figure 24 plots °3^ vs relative in-path slopes of from -30° 
to + 50° at 4 meters from the vehicle. In all cases, the cross 
path slope of the measured plane is 0° and there is at least 1 data 


point with a height h of 0 meters in each group of 4 points. 

.The graph only extends to -30° because the laser beam cannot 
"see" over a negative in-path slope greater than this value at 
4 meters distance, <5 r is varied from 1 cm to 5 cm with = 

1'. A data point spacing of .6 meters is assumed. Figure 25 
plots Cjg vis relative in-path slopes of from -6° to +50° at 20 
meters from the vehicle. All the other conditions are the same 
As in Figure 24 except that the data point spacing is 1.2 meters 
instead of 0.6 meters. Finally, Figure 26 shows graphs of 
vs relative cross-path slopes of from 0° to 50° for distances of 
4 meters and 20 meters from the vehicle. This time the in- path 
slopes are set to 0°. Once again 0/3 = <5© = 1' andO^t is varied 
from 2 to 10 cm. The data point spacing at 20 meters from the 
vehicle is 1.2 meters and the spacing at 4 meters from the vehicle 
is 0.6 meters. 

From Figures 24 - 26 it is noted that the higher the relative 
in-path or cross path slope, and the closer the data points are to 
the vehicle, the larger the error in gradient. All other factors 
considered equal, the error in gradient is larger for a given in- 
path slope than for a cross path slope of the same magnitude. It 
is obvious in Figure 24 that S'*, should be as close to 1 cm as 
possible in order that will be within 2° at 4 meters distance 
and high values of relative in-path slopes. 

This analysis leads to 3 very important conclusions. First, 
the very large error due to the motion of the vehicle is eliminated 
by utilizing a rapid scan laser. Secondly, if cf^ and ^ are re- 
duced below 1 arcminute, bulkier and heavier equipment would result 
In fact, with <3^ = <5$ of 1 arcminute, very little error is intro- 
duced into the calculations. Finally, most of the error is due to 
. The variance of the gradient should be as . low as possible, 
but if an upper bound of 2° is set, then Cr, must be as low as 1 cm. 

Projected Work Schedule 

Due to surface irregularity, data points will not fall at 
regular intervals as they would on a flat plane, but will be 
scattered about unevenly. A procedure must be developed to deter- 
mine how each group of 4 points is to be chosen prior to the 
modeling procedure. Another question that should be investigated 
is whether the modeling procedures should be modified for terrain 
at the 20-30 meter range so that the rover can model hills and 
craters as well as flat planes. More generally the problem is one 
of determining a complete picture of the terrain from an array of 
different sized and shaped planes, each of whose gradients have a 
different magnitude and direction. 
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d* Stereo and Non-Stereo Obstacle Detection Systems - W, Pfeifer 

Faculty Advisor: Prof. C. N. Shen 

The objective of this task is to evaluate obstacle detection 
systems based on alternative combinations of range and angle 
measurements. Both stereo and non-stereo concepts are being in- 
vestigated to determine the relation between measurement errors 
and the conclusions regarding the terrain feature. 

The concept is based on the assumption that a laser range- 
finder or its equivalent will be able to measure ranges and angles 
to terrain points on the Martian surface. The location of these 
points may be determined by several combinations of range and/or 
angular measurements. 

The Martian terrain has been modeled as a plane with discrete 
obstacles extending above and below it. Obstacles are defined as 
points on the surface with a height greater than .45 meters above 
or below the plane. Ref. 20. In addition, slopes between successive 
points must be considered because a change of height greater than .45 
meters spread out over a distance of several meters may be safely 
traversed by the vehicle. Both in-path and cross-path slopes must be 
determined since the vehicle has maximum limits of both pitch and roll 
from the local vertical. These limits are expected to be about twenty 
to twenty- five degrees. A method of fitting planes to data points and 
estimating the gradient of local surfaces is being developed under 
Task C. l.c. 

The coordinate system used in the investigations has the sensor 
mast as its vertical axis and a plane perpendicular to this axis at 
the surface. This plane may be referred to as a ,f vehicle base plane 1 '. 

Stereo Angles Method 

Preliminary investigation started with a two dimensional con- 
figuration which then can be generalized to three dimensions by the 
addition of one more measurement. A data point may be located by 
measuring the elevation angles from two knoton points on the vehicle 
to the point on the surface. This can be accomplished by using a 
transmitter (laser) at one location and a receiver at a second loca- 
tion on the vehicle. The transmitter and receiver are separated by 
a known distance L. As shown in Figure 27, by measuring the trans- 
mission and reception elevation angles, oc and /3 , and the known 
vehicle dimensions, the location (height above the vehicle plane and 
distance from the vehicle) of a point on the surface may be measured. 
Some type of scanning apparatus would aim the laser beam at a pre- 
determined pattern of points on the surface. 

The transmitter should be located at the lower point A as * 
shown in Figure 28 and the receiver at the upper point (B). With 
the transmitter at B, some data points will be lost since obstacles 
may intersect the reflected beam. 
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The values of h and b may be easily calculated once , /3 
and L and the remaining vehicle parameters are known from geometric 
considerations. 


L tan /3 
tan (3 - tan <*. 


( 1 ) 


tan - tan v ' 

The variances of h and b are computed from the variances of 
the angular measurements. It is believed that angular measurements 
with standard deviations of several arc seconds may be possible and 
that a standard deviation of one minute of arc is likely. Standard 
deviations of h and b have been calculated using the values of 6 arc 
seconds, 1 minute, and 10 minutes for the standard deviations of 
both oc and (3 . The value used for the standard deviation of L 
is .01 meters. 

The equations for variances are as follows: 
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Figure 29 shows the height standard deviation, &h , versus 
the distance from the vehicle, b, for h=0 and L=1.0 meters. Clearly 
the reduction of angular uncertainty from 10 minutes to 1 minute 
results in a significant improvement in accuracy. Also, further re- 
duction from 1 min. to 6 seconds results in slight improvement. Angle 
uncertainty should be reduced to 1 minute if possible, but further 
reduction is not very useful. Note that this stereo angles system 
determines heights very accurately over the entire 3 to 30 meter 
range when Cdf = 1 min. with being less than 4 cm. 


Varyinjg the value of the distance L between the transmitter and 
receiver also affects the accuracy of height determinations. Larger 
L gives better results, i.e. smaller . The percentage change in 
L results in about the same percentage change in The 1 meter 

value for L has been used in most calculations because the length L 
is limited by the vehicle design. Moreover, if point A moves too 
close to the surface the rangefinder loses some of its ability to see 
craters. 
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The distance measurement is not as accurate as the height. 

The distance standard deviation, <Tt> , versus distance, b, is shown 
in Figure 30 for h=0 and L=1.0 m. Sigma b is 1 meter for b greater 

than 16m and sigma angle = 10 min. With = 1 min., reasonable 

values of less than 12 cm are obtained f or at close range, i.e. 

less than 10m. 

It appears that the stereo angles method will be most useful 
at close range. At longer range, problems will be encountered be- 
cause of the inaccuracy in the determination of b. The distance 

determinations will be very important in the calculation of in-path 
and cross-path slopes. 

Stereo Rand Method 

Instead of measuring two angles, two ranges can be used to 
determine the coordinates of terrain points. To obtain two range 
measurements for a single terrain point a transmitter and receiver 
are located at one point on the vehicle and a separate receiver at 
a second point. Figure 31. Another possibility is to locate a trans- 
mitter between two receivers. 


The transmitter- receiver combination measures twice the dis- 
tance from the combination to the terrain point. The second receiver 
measures the sum of the distance from the transmitter to the terrain 
point and the distance from the point back to the second receiver. 

L, h and b are the same as in the stereo angles case. 

From the geometry in Figure 31: 


2 2 2 
; R = (3-h) + b 

2 2 2 
P = (3-h- L) + b 

Which can be solved for: 
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The actual quantities measured are defined as follows: 

M ■ 2R 
N * R + P 

If the range standard deviation is defined as , it 
follows that : 

Var M - 4 Var R 
* M .-2 

Var N = 4 Var R . 

<Sn = 2 <j R. 


(5) 

( 6 ) 

(7) 

( 8 ) 


(9) 

( 10 ) 


(ID 

( 12 ) 
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The variance of height and distance are then determined as 

2 

Var h - ^ lm ' L ^ Var L + ( 2 1 ) Var M + ( 2 "~ M ^ Var N (13) 

» - Uf [- 1 - ^ f L l H + N2m2 > } 2 V« L 

*uf 

♦(if 

The range standard deviation values used were 10, 6, 2, 1, 0.6, 
and 0.2 cm. Figure 32 shows height standard deviation versus dis- 
tance when h=0 and L=1.0 m. The height uncertainty is very large if 
range uncertainty is greater than 1 cm. To get <T h less than 10 cm 
at b = 7 meters CTj^ will have to be less than .5 cm. This type of 
range accuracy is not expected from present technology. 

The distance standard deviation, , versus distance is shown in 
Figure 33. The distance measurement has good accuracy when the 
range uncertainy isl centimeter or less. The distance uncertainty 
decreases slightly as b is increased from 7 to 28 meters. 

At present the stereo range system doesn T t appear to be useful 
because range accuracy is expected to be on the order of 5 to 10 
centimeters which leads to extremely large uncertainties in height 
determinations. However, if range standard deviation can be reduced 
to 0.5 centimeter or less this method would be acceptable. 

Non-Stereo (Range-Angle) Method 

This system, shown in Figure 34, uses a transmitter and receiver 
located at the top of the sensor mast. The rangefinder measures the 
range R to the terrain point at the elevation angle f3 to the point. 


From the geometry : 

h = 3 - R sin 3 (15) 

b = R cos & ; (16) 

Estimating the uncertainties in h and b using equations 19 
arid 20 leads to 

2 2 2 

Var h = sin £ Var R + R cos £ Var £ (17) 

Var b = cos (3 Var R + R sin Var (5 (18) 

Gh and crj, were calculated by assuming values for h, b, 


O'* and . Height standard deviation as a function of distance 
is shown in Figure 35 for range uncertainties of 10 era, 6 cm, and 
2 cm with angle uncertainties of 10 min., 1 min., and 6 seconds of 
arc. At close range, (less than 14 meters), range error seems to 


[ M - N 

l 2 

L 2 


(N M?N ) 
2L 2 


J 


+ ( 2N 3 -3N 2 M+NM 2 ) 


2L 


T~ 


] 


Var M 
2 


(14) 


Var N 


dominate. If angle uncertainty is less than some threshold value 
..between 1 and 10 minutes, height error actually decreases with 
distance at longer range. With <S* = 10 cm and = 10 minutes, 

is less than 9 centimeters over the entire 3 to 30 meter 
range. With angle accuracy of 1 min., height uncertainty is less 
than 4 cm even if <5^ = 10 cm. 

In the non-stereo system 5V was approximately equal to the 
range standard deviation over the 3 to 30 meter range. 

The non-stereo system is the most accurate method at ranges 
greater than 10 meters. If range error is 10 cm. and angle is 1 
min. then the non-stereo system is slightly less accurate than the 
stereo angles system for distances of less than 10 meters. 

The conclusion is that the non- stereo system is probably opti- 
mum because it has reasonable accuracy over the entire 3 to 30 meter 
range. 

Three Dimensional Systems 

The most logical way to extend the two dimensional systems to 
three dimensions is to measure an azimuthal angle 9 as shown in 
Figure 36. Coordinates a and b^ are .defined in the figure. 


a^ = R cos sin & = b sin 0 (19) 

b^ - R cos /6 cos 0 = b cos 0 (20) 

2 2 2 

Var a^ = sin © Var b + b cos © Var © (21) 

2 2 2 

Var b^ = cos 0 Var b + b sin B Var 0 (22) 


Equations 4, 14 or 18 may be substituted into equations 25 
and 26 to find Var a^ and Var b^ in terms of actual measured para- 
meters. The variances of a^ and b^ are expressed in terms of the 

variances of b and 0 . This means that the system which has the 

smallest variance of b will have the smallest variances of a^ and 
bi. If <*0 = 1 minute of arc, simple calculations show: 

if © = 0 ° cr El 1 cm 7 m < b < 28 m all C 

© = 12° 2.3 cm 7 i 4 b ^ 28 m OT) — 10 cm 

In general is approximately equal to . 

Effects of Heights Above or Below the Plane 

The stereo angle system is more accurate for positive heights 
than for negative ones as shown in Figure 37. For example, when 
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1 min. at a distance of 7 meters, °h = 2 cm if h = + 1 meter, 

= 3 cm if h = 0 m, and C* = 4 cm if h - - 1 m. This effect 
is certainly tolerable. Varying h has a negligible effect on 
in the stereo angles system. 

Non-stereo height accuracy is also better for positive heights 
than for negative ones. When angle standard deviation was 1 minute, 
the effect of. heights between plus 1 and minus 1 meter was similar 
to the effect on the stereo angles system and as range increased the 
effect decreased substantially as can be seen in Figure 38. Again 
this effect is tolerable. Varying h has a slight effect on <5^ at 
close range and almost no effect at ranges greater than 15 meters. 

In-Path and Cross-Path Slopes 

The in-path slope Sj and cross-path slope S e could be calcula- 
ted between successive data points. Ref. 21. . • 


S = 


A h h?-hl 

Ti = h£ ' - P j 


Var S T = 


Ab 

Var h 2 _ Var , 


2 1 


(23) 


(b 2 - h^ ^ Var + Var b ( 24) 

[_ ( Ab ) 4 l L 1 2J 


(JW- 

For cross path slopes the following substitution was included: 


Ab = [(a 2 - a L ) 2 + (b 2 - b x ) 2 ] h 

Var S = Var + Var h 2 ■ 

C ' TAbT 

O V 2 1 

+ [ ( h 2 " h i> \ \ ( a 2 " ap (Var a^Var a 2 ) + (b 2 -b ]L ) (Var bj+Var b 2 ) 

L if &)* J L tew (Ab ) 2 _ 

(25) 

Figure 39a shows the in-path slope standard deviations for 
the stereo angles system with 0^. equal to one minute of arc. The 
data point spacing, Ab , is 0.5 meter and 1.5 meters. The lines 
labeled S = 0° were determined from two points on the plane 
(hi « h 2 = 0). The lines labeled S = 25° were determined from point 
1 on the plane and point 2 a distance above the plane. 

h 2 = (Ab)(tan 25°) (26) 

The slope of a theoretical line between the two points is 25 degrees. 

When the two points are at the same height, the second term 
of the variance equation (24) is equal to zero which makes standard 
deviation of the slope smaller. The stereo angles system yields 
its lowest in-path slope accuracies at close range. 

Non-stereo in-path slope standard deviation is shown in Figure 39b 
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with Gr of 10 cm and of 1 min. Again two points at different 
heights give greater uncertainty in the in- path slope. With this 
system accuracy improves with increasing range. Non-stereo is 
more accurate than stereo angles for ranges greater than 10 meters. 

Cross-path slope uncertainties are nearly the same as the in- 
path slope uncertainties when the height difference is zero. Cross- 
path slopes do not appear to be very sensitive to the difference in 
height between the two points used in the calculation. 

Data point spacing should be increased as distance from the 
vehicle is increased. A larger number of points can be used in an 
averaging or plane fitting slope or gradient estimating scheme. 

Ref. 21. 

A rapid scan has been assumed so there is virtually no vehicle 
motion between the measurement of adjacent data points. The actual 
slope of the vehicle could be added afterward to the slope of the 
terrain measured to obtain the slope in relation to the local verti- 
cal. 

Conclusions 


Stereo angles may be a useful method, especially if range 
standard deviations are greater than 10 centimeters. The stereo 
angles method will be effective only at close range, i.e. within 
about 12 meters. The non-stereo method is almost as good at close 
range and is much better at ranges in excess of 12 meters. The 
stereo range method will not be accurate enough unless range un- 
certainty can be reduced drastically. The schedule for future 
work' during the next period is shown below. 

February March April May 

Investigation of 

data point spacing * 

Application of plane 
fitting to stereo 

systems — 

System evaluations 

and comparisons 

Final report -■ — — - 
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C.2. Path Selection System Simulation and Evaluation - S. Boheim 

Faculty Advisor: Prof. D. K. Frederick 

1b order to provide a means of evaluating proposed path selec- 
tion systems, the preparation of a comprehensive digital computer 
program was undertaken 18 months ago. The primary objective of the 
task has been to be able to simulate the operation of a vehicle on 
realistic three-dimensional terrains in the presence of random 
effects, with a variety of path selection systems. A secondary 
objective was been to establish quantitative criteria for evaluating 
a path selection system's performance, as simulated by the computer 
program referred to above. 

During the first academic year of work the general structure of 
the simulation was established and all of the essential elements were 
programmed, at least in an initial working version which was able to 
simulate one specific path selection system. The progress made in 
these two areas prior to July 1, 1972 has been summarized in Refer- 
ence 22. 

During the past six months, efforts have continued in the areas 
of simulation package improvements and path selection system design 
tradeoffs, and in compiling a cumulative technical report. Reference 
-23 and a user's guide for the simulation package. A brief discussion 
of these activities is presented below. 

A. Simulation Package Improvements 

Further attention has been directed towards the computer simu- 
lation package, as it was felt that a more powerful program would 
greatly aid in establishing improved evaluation criteria. The 
following additions and/or improvements were made to the simulation 
package : . . 

1. Vehicle Model - In the past, the vehicle has been 
modeled as a point source. The vehicle model has 
been updated to include the dimensions of the 
vehicle in some areas of the simulation package. 

Further work is necessary before this improvement 
is completed. 

2. Vehicle-Fixed Sensor - A sensor, with a fixed 
orientation with respect to the vehicle, has been 
simulated. This "vehicle-fixed" sensor is affected 
by the in-path and cross-path slopes of the terrain 
beneath the vehicle." 

3. Vertical-Fixed Sensor - In previous sensor simulators 
(see Ref erence 22) , the mid-range sensor which is as- 
sumed to have a range of 3 to 30 meters, was considered 
to be mounted on an ideal level platform directly above 


* The in- path slope is defined as the slope of the terrain measured 
in the direction of the vehicle's motion. The cross-path slope is 
measured in the direction perpendicular to the vehicle's motion. Both 
measurements are made by the vehicle. 


the vehicle. The "vertical-f ixed" sensor simula- 
tor has been improved so'that the effects of in- 
path and cross-path slope on the location of the 
sensor can be simulated. The vertical-fixed and 
vehicle- fixed mid-range sensors are compared in 
Figure 40. 

4. Range Measurement Simulation - Previous simulation 
methods (see Ref. 22) could calculate, with a 
resolution of only one meter, the length of the 
beam from the mid-range sensor to the beam impinge- 
ment point on the surface. Using this inaccurate 
range calculation as a first guess, an interval- 
halving algorithm has been utilized to obtain as 
fine a resolution as desired. For a typical de- . 
sired resolution of one centimeter, there was no 
noticeable increase in program execution time as 
compared to the time required to obtain one meter 
resolution. After further modifying the coding of 
sensor simulators, a saving of about 30% was 
realized in the execution time of these simulators. 

5. Short-Sighted Vehicle - The user has the option of 
not using a mid-range sensor on the vehicle. In 
addition, this feature allows simulation of a fail- 
ure of the mid-range sensor, if desired. 

6. Terrain Display Block - The terrain contour map has 
been expanded to display a larger area. The com- 
puter coding has been streamlined so that no appre- 
ciable execution time penalty for this size increase 
is realized. Two other maps have also been added. 

The first is a vehicle path map, which displays the 
path chosen by the vehicle, initial vehicle location, 
and the target location (see Figure 41). An appro- 
priate symbol is also used to indicate when the 
mid-range sensor is being used along the vehicle's 
path. The second map is an overlay of the vehicle 
path map onto the terrain contour map. 

7. Vehicle Dynamics Simulation - Noise may be added to 
the in-path and cross-path slopes of the terrain 
beneath the vehicle to simulate small terrain irre- 
gularities. Additional noise may be added to simu- 
late vehicle measurement error of these slopes. All 
noise addition is controlled by the user. The user 
may also limit the total simulation run time, thus 
preventing excessive costs when a path selection 
system becomes "lost” and wanders aimlessly (see 
Figure 41). 

8. Intermediate Output - To further speed execution 
time, the display of certain parameters of interest 
during a particular simulation run (such as vehicle 
location, heading, etc.), can be deleted by the 
user as desired. 
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B. Path Selection System Design Tradeoffs 

The performance of three path selection systems has been 
evaluated to form some preliminary conclusions regarding the 
tradeoffs in designing a path selection system. These conclus- 
ions (three of which are stated below) are considered in detail 
iii Reference 23. 

1. When using slope type terrain models, the 
obstacle thresholds (slopes are compared to 
these thresholds to determine the existence 
of an obstacle) should not be as large as the 
vehicle's maximum traversable slope. 

2. For slope type terrain models, increased use 
of the mid-range sensor allows better obstacle 
resolution, but may not be necessary if item 1 
is satisfied. 

3. Increasing the elevation angle of the sensor 
beam for slope type terrain models improves 
performance when in-path and cross-path slope 
effects are negligible. If these effects are 
not small, then the performance of the vehicle 
will be adversely affected. All of these 
situations are illustrated in Figure 42. 

C. Technical Report and User's Guide 

A technical report, Reference 23, has been prepared which 
summarizes all of the work done in this area of path selection 
system simulation and evaluation during the past 15 months. The 
development of the computer simulation package and the selection 
of evaluation criteria are discussed in detail. The performance 
of three different path selection systems is evaluated to deter- 
mine the effectiveness of the simulation package, and to form some 
preliminary conclusions regarding the tradeoffs involved in design- 
ing a path selection system. Some of these conclusions were presented 
in the previous section. Finally, suggestions for future work in this 
area are presented and are summarized in the next section. 

A user's guide has been prepared and will be available shortly. 
This text contains detailed information on how to use all of the 
available features in the simulation package, as well as the path 
selection systems presently implemented. Sections of the guide also 
specify guidelines for the user to follow when writing new path se- 
lection systems to be simulated and evaluated by the computer package. 
The derivation of all simulation formulas and methods used by the 
simulation package, as well as functional flowcharts, are included in 
the user's guide. 

Future Work 

1. Design Tradeoffs - Using the performance of three 

path selection systems, some preliminary conclusions 
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(A) Negligible In-Path and Cross-Path Slopes 



♦,><!), j no meaningful 
information available. 


(B) Significant In-Path and Cross-Path Slopes 


DESIGN CONSIDERATIONS FOR SENSOR BEAM ELEVATION ANGLES 

Figure 42 • 
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regarding the tradeoffs involved in designing a 
path selection system have been made. More exten- 
sive analysis is required, however, before more 
meaningful and reliable results may be obtained. 

Of particular interest is the comparison of the 
performance of the vehicle- fixed sensor with the 
vertical-fixed sensor, a study which will be con- 
ducted next semester and concluded by May. 

2. Evaluation Criteria - As path selection systems are 

evaluated, the desirability of these systems can be 
compared with the results of the simulation and 
evaluation package. A high degree of correlation 
is required between desirability and evaluation re- 
sults, certainly higher than is presently available. ’ '' 

Refinements and additions of evaluation criteria 

will be a continued effort that will become more 
effective as the number of evaluated path selection 
systems increase. 

3. Path Selection Algorithms - A large effort has been 
made to realistically simulate mid-range sensors, 
and, to a smaller degree, terrain modelers. On the 
other hand, only very simple path selection algorithms 
have been implemented. Work done next semester will 
involve the simulation and evaluation of more advanced 
path selection algorithms. 

Due to its modular arrangement, the program can be used to 
compare, by simulation, any number of path selection systems , regard- 
less of their specific features. Hence, it is anticipated that the 
program will be used extensively in the future to compare the per- 
formance attainable with sensors, terrain modelers, and path selec- 
tion' algorithms originating at JPL, Cornell, and other establishments 
in addition to Rensselaer. 

Task D. Chemical Analysis of Specimens 

One important phase of the initial missions to Mars is the search 
for organic matter and living organisms on the martian surface. The present 
concept for attaining this objective consists of subjecting samples of the 
atmosphere and surface material to certain chemical and biochemical reactions 
and thereafter analyzing the products produced, probably in a combination gas 
chromatograph/mass spectrometer (GC/MS). The gas chromatograph is proposed 
for separating complex mistures evolved from the experiments into small groups 
of similar chemical species. Chemical analysis of these groups would be 
accomplished in the mass-spectrometer. It is the objective of this task to 
provide engineering techniques and design criteria for designing such a system. 

The overall concept of the chemical analysis system is shown schema- 
tically in Figure 43. Most of the previous effort has involved the systems 
analysis of the gas chromatograph using simulation. Ref. 24, 25. This 
technique uses mathematical models, which incorporate fundamental parameters 
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evaluated from reported experiments, to explore various concepts and ito 
direct further experimental research. The objectives of the task have been 
extended to include preliminary studies of the entire chemical analysis 
system. System topics being considered include carrier gas generation and 
removal, and limitations imposed by and upon the mass spectrometer. 

The task problems are being attached by a four-member team, each 
of whom is pursuing a specific assignment: 

1. GC/MS system concepts 

a. Mass spectrometer system characteristics 

b. Carrier gas generation and removal 

2. Chromatographic system analysis 

a. Multicomponent chromatography 

b. Chromatograph model improvement 


D. 1. GC/MS System Concepts 

D. l.a. Mass Spectrometer System Characteristics - M. P. Badawy 

Faculty Advisor: Prof. P. K. Lashmet 

It is the function of the mass spectrometer to provide mass 
spectra for three different inputs (Figure 43) : effluents from the 

gas chromatograph; direct samples from the life-determining experi- 
ments; and samples of the atmosphere. This subtask has as its 
objective the characterization of the major properties of the mass 
spectrometer, the resulting relations to provide direction in de- 
signing an instrument suitable for this space mission. In particular, 
those characteristics influencing power, weight, and volume require- 
ments of the system are of primary concern in this study. 

Prior work provided a parametric study of the coupling between 
the magnetic field and electrostatic potential in determining the 
ion masses to be detected (Ref. 26). Instrument resolution and the 
nature of the ionizing process are currently being investigated. 


Resolution, or the ability of the instrument to differentiate 
among different masses at a given time, is estimated for the single 
focusing spectrometer as follows, Ref. 27, 28: 


M/AM = 



where 

s 1 , s 11 = slit widths at the source and detector ends, respectively. 


a^ = magnetic field radius of curvature. 

V = voltage used for deflection of mass M. 

£ V = voltage equivalent to range of masses emanating 

from the electrostatic field at different velocities. 
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^\M = dispersion in mass detected* 

Similar but more complicated results for the double focusing mass 
spectrometer have also been developed. These results emphasize the 
importance of the slit widths, the radius of curvature of the magnetic 
spectrum analyzer, and the energy spread of the instrument. Table 1 
^ summarizes the physical parameters of three mass spectrometers and 
presents their respective mass resolutions as estimated from the above 
equations. Higher mass resolution can be achieved with more complica- 
ted designs, e.g. the toroidal electric sector field design. Ref. 29. 
However, for the current preliminary studies, the characteristics of 
the more simple single or double focusing designs are satisfactory for; 
representing the basic instrument. 

Instrument sensitivity or the minimum concentration level detec- 
table is determined to a great extent by the ionization process. 

Ref. 34. The performance of the ionization chamber, based on an 
electron impact source, is characterized by the following equation. 

Ref. 30, 31. 

I = b Ncr s I 

e 

where 

I, I e - ** ion and electron currents, respectively 

b = ion extraction efficiency 

N = molecular density 

<3* = ionization cross-section 

s = effective electron path length 

Difficulty in determining the ionization cross-section limits the 
utility of the equation in the present study. Ionization cross-sections 
for single electron emission are available. Ref. 32. However, because 
knowledge of the components in the gas to be analyzed is uncertain, and 
determination of the ionization fragments is complicated, it appears 
that theoretical treatment of the ionization process is beyond the scope 
of this preliminary study. Experimental relations between the ion and 
electron currents will be sought to complete this phase of the task. 

Further work will be concerned with determining the effect of the 
major design variables upon power consumption of the mass spectrometer. 
Figure 44 shows schematically the relation between the various power 
consuming elements of the instrument. Current studies involve estimates 
of power associated with the filament of the ionization chamber and with 
maintaining the accelerating potential. After these studies are com- 
pleted, the vacuum pumping system and the magnet power supply will be 
considered. 

D. l.b. Carrier Gas Generation and Removal - C. W. Jarva 

Faculty Advisor: Prof, L. K. LAshmet 

Carrier gas, which is relatively inert with respect to the samples 


84 


'll 


p 

<D 

CO 

0 

■fi 

O 

O 

> 


*6 

o\ 


o 

VO 


oo 


* 

o 

+1 


8 • 6 


-=t ' 

ir\ 

'iR 

. 

CVJ 

ir\ 

CO 

cvJ 

• 

to 

• 

H 

€ 

ir\ 

rH 

CVJ 

• 

• 

rH 

CVJ 

oo 

O 

O 

H 


. 


d 

o 




co 

fH 

0 

P 

0 

I 

*P 

o 

a> 

CO 

CO 

CO 

a5 


SI 


*d 

§ 

S 

1 

w 


1 

1 

8 cm 

0.02* 

g 

f; 

g 

Kr 

g 

■ g 

CO 

oo 

O 

oo 

VO 



• 

o 

+ 1 

CVJ 

• 

-3“ 

* 

o 

• 

o 

CVJ 


iH 

oo 


W 


H 


a? 

P 

0 

£ 

CO 

p 

o 

CO 

p 

<u 

-p 

<L) 


l-H 

Ph 

^>* 


8 


o 

o 

VO 

LTV 

• 

0 

CVJ 

• 

« 

VO 

o '■ 

CO 

o 

o 

• 



OY 

o\ 

H 

+1 

oo* 




& & 


o 

O 


ITN 

o 

o 


VO 

H 


OV 

o 

OO 

H 


aj 

Ph 

O 

<P 

•H 

O 

0 

& 


<D 


'SC 


0 


a) 

•rH 
<H 

O 

•rH 

P _ 

4? •A 

w d 

8 05 
P 
O 
0 

d 


0 

"A 

s 


o 

•H 

P 

a5 

P 

to 


5 


aJ' 


S 


0 


o5 


*d 

O hO 

o 

*d 


rH 

0 a 

•H P 

rH 

o 

0 

H 0 

CO £ 

0 

•H 

•H 

0 H 

P 0 

•H 

P 

<H 

o *d 

0 Ti 

<H 

d 

Pi O 

CO 

P 

o 

P H 

co »rl 

a d 

co 

•H 

P 

0 

0 *H 


•H *H 
P -d 

P 

0 

d 

O <P 

0 

0 ctf 

P 

P 

co o 

d P • 

o 

ho 

d 

co o 

hb 

0 

3- 

o 

co 

d 

£5 

1 

a 


CO 


*d 

"d 

•H 


o3 

fn 


I 


P 

•H 

rH 

CO 

0 

O 


P 


CO 


g 


p 

•H 

H 

co 


Ph 

CO 


£ 


g 

<3 

s' 


g. 

•H 

P 

3 

O 

CO 

0 

u 

co 

CO 


'd- 

0 


co 

3 


STORAGE 
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being analyzed and which must be in supply for an appreciable part 
of the mission, is required to propel minute samples through the 
column of the gas chromatograph, Ref. 33. On arriving at the mass 
spectrometer, which operates at a low pressure ( < 10“^ torr) , the 
major portion of this carrier gas must be separated from the sample 
and must be removed from the system. Thus the systems for genera- 
ting, separating, and removing the carrier gas are important to 
overall system design and operation. It is the objective of this 
subtask to generate the basic characteristics and limitations of 
several promising concepts to aid in the design of the final overall 
system. 

A study was undertaken to determine, as a function of the initial 
charge pressure, the physical characteristics of a spherical, high- 
pressure gas tank for storing the amount of carrier gas necessary for 
experimentation. Because of diffusion of the carrier gas, assumed to 
be hydrogen, through the vessel wall and its leakage through welds, 
valves, etc., Ref. 34, 35, the gas storage capacity was greater than 
the requirement for the experimentation. Diffusion rates were deter- 
mined from reported data, Ref. 36. To estimate leakage through welds 
and valves, it was assumed the vessel would be subjected to helium 
leak detection when charged to a practical pressure taken as 34 atmos- 
pheres in this case. The minimum leak detectable by modern leak de- 
tectors, about 10" 6 atm, cnrVsec, Ref. 35, was extrapolated to account 
for hydrogen as the gas and for the actual charge pressure. The vessel 
wall thickness was estimated from standard equations, Ref. 37. 

Estimates of vessel size and weights as affected by charge pres- 
sure and operating time are shown in Figures 45 and 46. These estimates 
are based on the following assumptions: 

1. Construction - fiber reinforced stainless steel 
containing 1/3 metal by volume, Ref. 38. 

8 2 

2. Allowable working stress - 4.5 x 10 N/m 
(66,000 psi), Ref. 38. 

3 

3. Gas flow rate for experiments - 2 std cm /min, Ref. 39. 

4. Flight time to Mars - six months. 

5. Maximum temperature - 300°K. 

Although. diffusion and leakage were considered in these studies, these 
processes contributed only about 1% to the vessel volume. The vessel 
weights, shown in Fig. 46 to be approximately independent of charge 
pressure, are for the shell only. The pressure regulator and its 
mounting boss will contribute an additional one kgm approximately. 

It was a major objective of this study to examine the merits of 
the gas storage vessel especially relative to a water electrolyzing 
cell proposed earlier. Ref. 40. It is estimated that the gas storage 
tank will weigh approximately 60% more than a comparable electrolyzing 
cell system but will occupy about the same space depending upon the 
charge pressure. It thus appears that if reliable fiber- reinforced 



Figure 45. Effect of Cha 
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vessels having about four times, the allowable working stresses of 
stainless steel, can be fabricated, the concept of the high-pressure 
storage vessel has promise, especially since the system is simple. 

Except for possible minor improvements in the calculational 
procedures, this preliminary study has been completed. 

A carrier gas removal system must be considered because the 
operating pressure of the mass spectrometer is about 10“^ times the 
pressure expected in the chromatograph. Ref. 41. The palladium- silver 
alloy/hydrogen gas separator has been proposed. Ref. 39. In such a 
system, hydrogen flowing in a palladium tube diffuses rapidly through 
the heated tubing wall in comparison to other gaseous materials. This 
subtask is concerned with developing technical procedures suitable for 
the design of such a unit as well as with assessing its merits. 

At present the steady state behavior of the system is being con- 
sidered using momentum balance and material conservation equations. 
Hydrogen diffusion rate through a palladium wall is proportional to 
the difference of the square roots of the hydrogen gas pressures at 
each face of the wall because of dissociation during diffusion, Ref. 

42. A diffusivity particular to the system is also involved, Ref. 42. 
From this diffusion rate equation, a first order, nonlinear, differen- 
tial equation has been developed which relates the rate of material 
flowing at any point in the tube to the rate the gas is diffusing out 
of the system. Combined with a second order nonlinear momentum balance 
equation, a describing set of equations for this system was formed 
which requires numerical solution. This analysis assumes laminar flow, 
pure hydrogen flowing through the palladium tube (no presence of sample), 
and radial flow of the diffusing gas through the palladium walls. Equa- 
tions accounting for the presence of sample will be developed after the 
characteristics of this simple system have been determined. 

In the present studies, the hydrogen pressure outside the palla- 
dium tube is assumed to be approximately zero. Pumping means such as 
combining the diffused hydrogen with oxygen on the tubing wall. Ref. 29, 
is a topic of future work. Other topics worthy of consideration . if 
time permits include estimations of mechanical breakdown due to tempera- 
ture and pressure cycling and chemical breakdowns due to contaminants. 



D.2.a. 


Multicomponent Chromatography - A. J. Meisch 
Faculty Advisor: Prof. P. K. Lashmet 


Linear mathematical models of the chromatograph have been developed 
for single chemical components and have been solved analytically for an 
impulse input, Ref. 25, 43, 44. Other work. Ref. 24, which used data 
from a chromatographic test facility has shown that convolution of the 
actual input data with the theoretical normalized impulse response gives 
a reasonable representation of the actual chromatogram for a single 
component. 
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It is the objective of this subtask to experimentally evaluate 
the mathematical model for multicomponent chemical systems using 
superposition of a single component data. Specifically the equili- 
brium adsorption model, Ref. 44, is being considered: 

y (e) = i \l/3p e /4 rre^exp [_ - (P e /4/je) (e -ft ) 2 

where 

(E = 1 + (1/m R o ) 

0 ~ dimensionless time = vt/L 

v = carrier gas velocity 

t = time 

L = column length 

Pe = the Peclet number, a dimensionless measure of 

sample diffusion in the carrier gas 

ihRq - a thermodynamic parameter, specific to the chemical 
system considered 

Y = composition response to a unit impulse. 

This model assumes the gas at each point to be in equilibrium with 
the solid adsorbent; i.e., the column is very long. The Peclet number 
is a function of the fluid mechanics and the physical properties of 
the system and is predictable, Ref. 25. The parameter hiRq is specific 
to the system and is determined from the system data using a curve 
fitting technique, Ref. 27. 

Extension of the model to multicomponent chemical systems using 
superposition of pure component data was proposed. Experimental work 
reported earlier, Ref. 25 on the binary system n-pentane/n-heptane on 
a microporous adsorbent Chromosorb 102 column (see Table 2) showed 
superposition was only a first order approximation and large deviations 
were noted. This subtask has as its objective the verification of the 
previously observed discrepancies and the generation of additional 
multicomponent chromatographic data. 

Chromatographic data for the system n-pentane/n-heptane on the 
Chromosorb . 102 column have been obtained at 150 and 200°C and for 
mixtures containing 0, 1, 25, 50, 75, 99, and 100% n-pentane by weight. 
As shown in Figure 47, the peak times for the multicomponent system 
lagged peak times predicted from superposition of the pure component 
data. This time lag occurred for both compounds for all mixtures con- 
sidered as shown in Table 3. The time lag increased as the component 
composition became more dilute. 

The peak time is primarily related to the parameter mRg, Ref. 44. 
As the peak time varies with sample composition, it is possible that 
mRo is composition dependent; ;i.£. , although the component compositions 


LENGTH (CM) 

0 D 

1 D 
MESH 

PARTICLE SIZE 
COMPOSITION 

TEMPERATURE RANGE 
APPLICATION 


TABLE 2 

COLUMN CHARACTERISTICS 

CHROMOSORB 102 

100 CM 
i/8 IN 
0.22 CM 

6o/8o 

0.025/0.0177 CM 

MICROPOROUS 
STYRENE-DIVINYL 
BENZENE POLYMERS 

ROOM -250°C 


SEPARATION OF LOW 
MOLECULAR WEIGHT, 
VERY POLAR SUB- 
STANCES SUCH AS 
WATER, ALCOHOLS, 
AMINES, ACIDS, 
ESTERS, KETONES, 
ETHERS, ALDEHYDES 
AND TO A LIMITED 
EXTENT LIGHT GASES. 


POLYETHYLENE GLYCOL 


100 CM 
1/8 IN 
0.22 CM 
60/80 

0.025/0.0177 CM 

POLYETHYLENE GLYCOL 
(CARBOWAX 1500) 

20$ BY WEIGHT 

ROOM-22 5°C . 


SEPARATION OF HIGH 
BOILING POLAR COM- 
POUNDS. GENERAL 
PURPOSE COLUMN FOR 
POLAR SAMPLES.. 


PLOT OF PLKF CCFPCAENT OAT A VS. ACTUAL l.‘ A T A 
^ACTUAL PENTANE DATA, FACTUAL HEPTANE DATA, 2 = PLKfc PtMANfc. LATA, l*Pl,^E FtPTANE CATA 



TABLE 3 

PENTANE - .HEPTANE SYSTEM ON A CKROMASOKB 102 COLUMN 
HELIUM FLOW HATE OF 29.4 CC/MIN AT 150 DEG C 



ACTUAL 

ACTUAL 

COMPOSITION 

PENTANE 

HEPTANE 

(WGT <f>) 

PEAK * 

PEAK * 

100 H 

• 

285.5 

99 H - IP 

95.. o 

294.0 

75 H - 25 p 

88.5 

307.5 

50 h - 59 P 

86.0 

327.5 

25 H - 75 p 

84.0 

364.0 

1 H - 99 p 

81.5 

432.5 

100 P 

80.5 

•a 


* TIMES ARE IN SECONDS 
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in the helium carrier gas are small (^10 ^ mole fraction), there 
are possible interaction effects at the adsorbent surface. 

Future work on the n-pentaiie/n-heptane system will include 
data generation at 175° C on the Chromosorb 102 column and at 75° 
to 125° C on a nohporous Carbowax 1500 column described in Table 2. 

The chromatographic test facility described earlier, Ref. 45, 
is being modified to improve data handling. Automatic data proces- 
sing as outlined in Figure 48 is being implemented to improve con- 
sistency in the data and to reduce the manual effort in handling. 
Oscillographic strip charting will be continued to verify the output 
from the more automated procedure. 

Studies have also begun on the basic characteristics of the 
adsorbents used in the chromatographic columns. This work includes 
experimental determination of pore size distribution in the adsorbent 
particles and void fractions in the columns. These data are required 
for accurate prediction by the improved mathematical models being 
developed (see Task D.2.b.). 


D.2.b. Chromatograph Model Improvement - P. T. Woodrow 

Faculty Advisor: Prof. P. K. Lashmet 

Prior to July, 1972 the mathematical model of the chromatographic 
column was based on a packing material which either was non-porous or 
which allowed pore diffusion and subsequent adsorption to occur at 
such a high rate that the resultant system dynamics were not signifi- 
cantly affected. Thus, in modeling the system, only interp article 
and adsorbed phase mass balances were considered. The applicable partial 
differential equations, assumptions, initial conditions, and boundary 
conditions have been presented elsewhere. Ref. 25. These equations 
showed the predicted chromatogram to be a function of three variables: 

jiiRq = a thermodynamic parameter, specific to the chemical 
system considered. 

N nr = the number of transfer units, a dimensionless 

J measure of diffusion rates to the adsorbent surface. 

Pe = the Peclet number, a dimensionless measure of 
sample diffusion in the carrier gas. 

Solutions to this system of equations have been reported, and 
cases of equilibrium adsorption (y = y* ; N^qq~> oo ) and non- 
equilibrium adsorption (y = y*; N t0 Q is finite) have been compared 
with experimental data, Ref. 25. It was concluded that for the ranges 
of parameters encountered in typical chromatographs, consideration of 
a finite rate of mass transfer (NtOG = finite) yields minor, at best, 
improvement in output chromatogram predictions. In addition, the 
predictions of peak height and spreading were inadequate for some 
experimental conditions and chemical systems. It appeared that another 
mechanism not considered in the model formulation was responsible for 
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the discrepancies between prediction and experiment* As the adsor- 
bent (Chromosorb 102) was porous, Table 2, the mechanism of intra- 
particle diffusion warranted investigation. 

Consideration of the intraparticle diffusion process as well as 
the rate of surface adsorption yielded a new system of equations 
shown in Fig. 49 and having the following parameters and variables 
not defined, previously: 

thermodynamic equilibrium constant. 

the number of reactor units, a dimensionless 
measure of the rate of adsorption. 

intraparticle Peclet number, a dimensionless 
measure of diffusion rates within the particle. 

interparticle Peclet number, a dimensionless 
measure of diffusion rates in the carrier gas. 

ratio of column length to particle radius. 

moles of fluid in particle per mole of adsorption 
sitqs. 

intraparticle space variable, dimensionless. 

adsorbed phase composition, dimensionless. 

interparticle gas phase composition, dimensionless. 

intraparticle gas phase composition, dimensionless. 

equilibrium intraparticle gas phase composition, 
dimensionless. 

z = longitudinal position in column, dimensionless. 

As an aid to envisioning the various transport processes incor- 
porated in this new chromatograph system formulation. Figure 50 is 
included. In the interparticle region, material is transported by 
convection (i.e. , by carrier gas flow) and by turbulent and molecular 
diffusion in the axial direction, which is characterized hy the inter- 
particle Peclet group. Peg. The transport from the interparticle 
region to the intraparticle region is characterized by the dimension- 
less mass transfer coefficient, N t0G . Transport by diffusion within 
the particle - the intrap article region - is characterized by the intra- 
particle Peclet group, Pe^. The rate of adsorption from the intraparticle 
gas phase to the adsorbed phase within the particle is characterized by 
Nru> fc h e number of reactor units. 

The transfer function of the system of Figure 49 is given in 
Figure 51. Before attempting an inversion of the equations of Figure 51, 
a moment analysis of the model hereafter, designated as the Inter- 
Intraparticle Adsorption Model was undertaken to determine gross behavior. 


m 


N, 


RU 

Pe. 


Pe E = 

L/R = 
R_ 


y 

y. 


y? 


interparticle phase mass balance: 



Figure 49. System Equations for Inter 
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Figure 51. Transfer Function of Chromatograph 
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This technique has been described earlier. Ref. 44. 

Moments derived from the equations of Figure 51 provide 
characteristics of an impulse response. Because the inputs to the 
system are pulses of finite duration, equations relating the moments 
of the output chromatogram with the moments of the impulse response 
and input were derived : . 


\ \ 

Mi y = + ® . 

A >2 = >2 + >2 

Y X y 

= >3 X + f 3y 


where 


V'St 

\A 


areas under the output and input curves, 
respectively. 

first moments about the origin of the 
output and input curves respectively. 


6 



first moment about the origin of the 
impulse response. 

second moments about the mean of the output 
and input curves respectively. 



second moment about the mean of the impulse 
response. 



third moments about the mean of the output 
and input curves respectively. 



third moment about the mean of the impulse 
response. 


The three moments 9, , and p were derived from the equations 

z y __ _ 

of Figure 51, whereas the input moments Jw~ and Icl 

' ly, 2 X , ' 3y were com- 


puted directly from experimental data. A priori prediction of XL/R), 
P ®A> ^ e E> an< * ^toG was P oss: *-ble from the experimental operating con- 
ditions. The parameters mR^ and must be determined by curve- 

fitting procedures. 


To determine qualitatively the capability of the new model of 
Figure 49 fore representing more adequately the chromatographic data,, 
a parametric study of the moments was undertaken. Table 4 compares 
moments obtained for the ethylene system. Ref. 25, with predictions 
from the interparticle equilibrium adsorption model. Ref. 25, 44, and 
the new model. The new model more closely predicts the mean, yhq, 
although differences between the two models and the data are small. 


lor 


TABLE 4 

MOMENT ANALYSIS AND PARAMETRIC 
STUDY - ETHYLENE 50 °C. 




observed 



(i) 

predicted 



( 2 ) 

predicted 


0.194 26.475 25.986 23.719 




observed 


100 

200 

400 

800 

1600 

32OO 

6400 

12800 

25600 


7.024 


1 




predicted 


-( 2 ) 

predicted 


13.283 

6.973 

3.817 

2.240 

1.451 

I.O56 

0.859 

O.76O 

0.711 


0,388 


i 



w 3 

observed 


3 

predicted 


“( 2 ) 


^predicted 


100 

200 

4oo 

800 

1600 

3200 

64oo 

12800 

25600 


19.623 

-L 


13.049 

3.519 

1.058 

0.403 

0.219 

0.163 

o.i44 

0.137 

0.134 


0.191 


± 


Pe. 


9744 


N 


toG 


79750 


(l/r) /P e = 436.2 

A 


1) Inter-Intraparticle Adsorption Model 

2) Interparticle Equilibrium Adsorption Model 
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More importantly, the new model is capable of better predicting the 
dispersion, as represented by the second moment about the mean,^, 
and more closely approximates the skew as represented by the third 
moment about the Values of necessary for good 

agreement between the moments are in reasonable accord with other 
independent data, Ref. 46. It thus appears that the proposed model 
will yield more adequate chromatogram predictions. 

The solution to the equations of Figure 49 is presently being 
considered. It appears that inversion of the transfer function of 
Figure 51 is not feasible, and an analytical solution is probably not 
realizable. The simplified situation of negligible mass transfer 
resistance (N t oG->c>o) and equilibrium adsorption was also 

investigated. It appears the resulting transfer function, while con- 
siderably less complicated, is not easily inverted. Future activities, 
therefore, will consider direct numerical solutions to the system 
equations. . 
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